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ABSTRACT: We review the current status of determinations of the b-quark mass, mb- We 
describe the theoretical tools required for determining nib, with particular emphasis on effective 
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Cu field theories both in the continuum and on the lattice. We present several definitions of mb 

and highlight their advantages and disadvantages. Finally, we discuss the determinations of rub 
from bh systems, 6-flavored hadrons, and high-energy processes, with careful attention to the 
corresponding theoretical uncertainties. 
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1 INTRODUCTION 

The mass of the bottom quark, nif,, is a fundamental parameter of the standard 
model. As such, it is an important quantity to measure, both for its own sake and 
as an input into the determinations of other parameters. However, because of con- 
finement, nib (like any quark mass) is difficult to determine experimentally. Free 
quarks do not exist; hence, rrn, must be inferred from experimental measurements 
of hadron masses or other hadronic properties that depend on it. Furthermore, 
like any parameter in the Lagrangian of Quantum Chromodynamics (QCD), rrih 
is a renormalized quantity, and therefore scheme- and scale-dependent. Although 
any renormalization scheme is possible in principle, some schemes are more con- 
venient for a given purpose than others. 

In this review, we discuss the current status of rrih. In Section |2|, we briefly 
introduce the main theoretical tools in current use: effective field theory (EFT), 
lattice field theory, and the operator product expansion (OPE). In Section ^ we 
discuss the relative advantages of several popular quark-mass definitions, and in 
Section ^ we discuss the various determinations of nii,. 

1.1 The Importance of nib 

In the standard model, quark masses arise from the coupling of quarks to the 
Higgs field, which acquires a nonzero vacuum expectation value through sponta- 
neous symmetry breaking. These couplings are all free parameters in the standard 
model, so a precise determination of rrib is interesting both in its own right as a 
fundamental parameter in the standard model and as a constraint on models of 
flavor beyond the standard model (see, e.g.. Reference ([l|)). 

From a more phenomenological perspective, a precise value of rrib is an impor- 
tant ingredient in the current experimental program of precision flavor physics 
at the B factories (and elsewhere) (^). Studying B decays and mixing allows 
the elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix 
to be over constrained, and hence provides a sensitive test of the flavor sector of 
the standard model. Because the predictions for many standard-model processes 
are very sensitive to mf,, the uncertainty on rrib feeds into the uncertainties of 
other parameters, limiting the precision at which the standard model may be 
constrained. An important example is the rate for the inclusive semileptonic 
b — > uiv decay used to determine the CKM matrix element \Vub\-, which is pro- 
portional to m\. The determination of \Vcb\ from inclusive h — > civ decay is less 
sensitive to mb than \ Vub\ is, but rub is still an important source of uncertainty. 

Using heavy-quark symmetry, a determination of the 6-quark mass also yields 
a value of the charm-quark mass. An accurate determination of the charm-quark 
mass, mc, is also required for precision tests of the standard model. For example, 
the rare decay — > tt'^vu is sensitive to the CKM element \Vtd\ and depends 
on rric through virtual charm loops (^). 

Finally, the different determinations of mb use most of the tools that have 
been developed to deal with the physics of hadrons: heavy-quark effective the- 
ory, lattice QCD, and operator product expansions (OPEs) applied to inclusive 
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decays. These are applied to observables as disparate as the masses and widths 
of low- lying T states, the masses of the B and Bg mesons, and the moments 
of inclusive i?-meson decays. Different determinations have completely different 
sources of theoretical and experimental uncertainty, so the consistency between 
these different determinations is an important test of our theoretical tools. 



1.2 Model Dependence and Theoretical Errors 

The earliest values quoted for the 6-quark mass were determined by fitting the 
observed spectra of hadrons to phenomenological quark-antiquark potentials (^, 
^). A simple constituent quark model in which a ground-state 6-flavored hadron's 
mass is given by the constituent quark masses and a spin-spin interaction 

M = J2m, + aY^^l^ (1) 

I i<j J 

is remarkably successful at reproducing the masses of the ground-state hadrons 
and it leads to a constituent mass mb ^ 5 GeV. Similarly, a simple linear- 
plus-Coulomb-potential model @, 

V(r) = -- + ^ (2) 

(where k is some effective strength of the potential), leads to m^, ~ 5.17 GeV from 
the measured T{2S) — T(IS') splitting. 

Because neither Equation [l| nor Equation |2| is derived from QCD, these de- 
terminations of nib carry serious disadvantages. The constituent 6-quark masses 
in these equations are defined only in the context of a quark model and have 
no well-defined relationship to the quark mass in the Lagrangian. Hence, the 
resulting determination is model-dependent, and it is difficult to assign a sen- 
sible theoretical error to the result. Much refinement of the simple potential 
model in Equation ^ is possible (|5|), but such refinements do not parametrically 
bring the potential model closer to QCD and do not systematically improve the 
determination of nib- 

Control over the theoretical errors is an important issue. Ultimately a precise 
determination of rrib is interesting because, along with other precision measure- 
ments, it may teach us something about new physics. If a discrepancy arises 
between precision data and theory, we can claim that the data point toward new 
physics only if the theoretical errors are truly under control. Therefore, this 
review concentrates on model- independent determinations of nib- We define a 
model-independent result as one that becomes exact in some limit of the the- 
ory, so that the theoretical error is quantifiable (determined by some parameter 
that describes how close the real world is to this limit) and in principle may be 
systematically improved upon. 

Although the size of theoretical errors from model-independent results is para- 
metrically known, one must often still resort to models (if only dimensional anal- 
ysis) to estimate the size of the uncertainty. There are a number of sources of 
theoretical uncertainties. The continuum calculations discussed below contain 
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both perturbative and nonperturbative uncertainties, proportional to powers of 
as{mb) and AqcD/n^b, respectively. In lattice QCD calculations, theoretical er- 
rors arise from the approximations used in the numerical simulations: discretiza- 
tion effects due to the finite lattice spacing, the incomplete inclusion of sea-quark 
effects (quenched approximation), and heavy or light quark-mass extrapolations. 
However, the uncertainties associated with these effects can, in principle, be deter- 
mined from within the lattice calculation by varying the parameters and studying 
their effect on physical quantities. Perturbation theory, if used in both lattice 
and continuum calculations, gives rise to theoretical errors due to the truncation 
of the perturbative series. Because perturbative uncertainties rely on an assump- 
tion about the size of the uncomputed higher-order terms, they depend on the 
quality of the perturbative expansion. 

Theoretical errors must therefore be treated with caution, akin to systematic 
errors for experimental results; they represent the theorist's best guess at the 
expected size of uncomputed terms, and it is often difficult to quantify the un- 
certainty of the error. Finally, when quoting determinations of m;, with different 
sources of theoretical error in this paper, for brevity, we combine these in quadra- 
ture to quote a single theoretical uncertainty. 



1.3 The Heavy- Quark Limit 

It is very useful to consider the limit of QCD in which the 6-quark mass is much 
larger than the QCD scale, 

0. (3) 



Aqcd 



rrib 

In this limit, the physics that occurs at a distance scale l/m;, may be cleanly 
separated from the physics of confinement, which occurs at the much larger dis- 
tance scale 1/Aqcd- This has several benefits. In some cases, the hadronic 
physics in a process of interest may be related through approximate symme- 
tries to the hadronic physics in another easily measured process. In others, the 
hadronic physics is irrelevant in the heavy-quark limit, and the process is purely 
determined by short-distance physics. The latter situation is useful for our pur- 
poses because short-distance physics probes the b quark properties, whereas long- 
distance physics is sensitive to the b hadron properties. Hence, physical quantities 
that are determined by short-distance physics provide a sensitive probe of the b- 
quark mass. A familiar example, which we discuss in more detail below, is the 
i?-meson semileptonic width. In the limit of Equation |^, this is a purely short- 
distance process, and, in this limit, the total semileptonic width of a meson is 
the same as that of a free quark. 

The heavy-quark limit is not a bad approximation for real b quarks. Because 
^QCD/nT-b ~ 0.1, we would naively expect corrections to this limit to be of order 
10%. Nevertheless, for precision physics, the corrections of order and 
(Aqcd/'^^)^ (if not higher) should be taken into account. This becomes a com- 
plicated problem. Not only must all finite-mass effects be accounted for order 
by order, but virtual loops probe both long- and short-distance scales: Radiative 
corrections therefore introduce additional corrections to the heavy-quark limit. 
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which are suppressed by powers of as{m{j). Furthermore, other energy scales are 
relevant to b decays, particularly the charm-quark mass nic- For the T system, 
both the three-momentum nihV and the kinetic energy mf,v'^ of the b quark enter 
the dynamics. The simplest way to keep all of this under control is by means of 
effective field theory. 

In the next section, we discuss several EFT approaches that are useful for 
determinations of m^. 



2 EFFECTIVE FIELD THEORIES, OFF AND ON THE LAT- 
TICE 

Effective field theory (EFT) is a general tool for dealing with multiscale prob- 
lems by separating the contributions from the different scales (for several lucid 
reviews, see Reference (0)). Interactions due to short-distance physics can, in 
general, be described by local operators. The coefficients of these operators de- 
pend on the short-distance physics and are usually calculable in perturbation 
theory by matching calculations of physical (on-shell) quantities in the EFT to 
their counterparts in the underlying theory. 

In our case, the EFT approach allows us to separate the short-distance physics 
associated with the the 6-quark mass from the long-distance QCD dynamics. 
The idea is that the leading-order Lagrangian corresponds to the heavy-mass 
limit; corrections to this limit are incorporated by nonrenormalizable operators, 
suppressed by powers of l/rufj. The theory can be renormalized order by order 
in 1/mf,. It is important to note that despite the label "effective," any quantity 
calculable in the full theory is calculable (usually much more simply) in the 
effective theory. 

Because we apply the EFT to the b quark only, the Lagrangian takes the form 

^cS = -^heavy + -flight, (4) 

where /^heavy is the effective Lagrangian for the b quark, and £iight is the usual 
QCD Lagrangian for the gluons and Uf light quarks, 

Aight = -^Tr G^.G'^^ + J2 QiUP - mi)qi. (5) 

i 

In the following two subsections, we discuss two different forms of £heavy) which 
correspond to the heavy-quark effective theory (HQET) and nonrelativistic QCD 
(NRQCD) treatments of the b quark. 

The EFT approach is not only important in continuum calculations but also 
essential in lattice calculations. It allows us to separate the short-distance effects 
of the lattice discretization from the long-distance nonperturbative dynamics and 
hence to use field-theoretical methods to analyze and correct discretization effects. 



We discuss strategies for dealing with heavy quarks in lattice QCD in Section 2.3. 
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2.1 HQET 

HQET was developed primarily in the late 1980s and early 1990s (^J^, IC , 1 1 , 1^, 13 ) 
as a tool to systematize the simplifications of the heavy-quark limit and to sim- 
plify calculations of both perturbative and power corrections to this limit. Many 



excellent reviews of HQET (see, e.g., Reference (14)) provide a more complete 
discussion. 

The effective theory is constructed by splitting the momentum of a heavy quark 
into a "large" piece, which scales like the heavy-quark mass mg, and a "small" 
piece (the "residual momentum"), which scales like Aqcd: 

= rngv" + k". (6) 

Here mg may be any mass parameter that differs from the meson mass by a 
term of order Aqcd- Typically, it is chosen to be the pole mass, discussed in 
more detail in Section |3.1| , but other choices are possible [and equivalent (|l5|)]. 
The two terms in Equation ^ are distinguishable because, in the limit uiq — > oo, 
interactions with the light degrees of freedom cannot change but only k^. 
The four- velocity v^^ is therefore a conserved quantity in the EFT (p!3[), and may 
be used to label heavy-quark states. The HQET Lagrangian is defined as an 
expansion in powers of k^/rriQ ~ Aqcd/"^q, 

-^heavy = Cq + Ci + C2 + ■ ■ ■ , (7) 

where Cn is of order l/rriQ. The first few terms are 

£0 = hyiD ■ vhv (8) 

and 

A = ^K{iDfK + ^^Ka^^G^pK, (9) 
2mQ 4mQ 

where hy is a heavy-quark field with four- velocity v, Ck = 1 (to all orders in 
perturbation theory (|l6|)), and CmifJ-) = 1 + 0[as{mQ)] is known to two- loop 
order (p^. The leading term. Equation!^, is spin- and fiavor-independent. Hence, 



at leading order, the interactions of heavy quarks with light degrees of freedom 
have a spin-fiavor symmetry, which dramatically simplifies the study of exclusive 
decays, as well as hadron spectroscopy, in the heavy-quark limit (^, p, IC, Tl]). 



The two terms of order l/mq correspond to the heavy-quark kinetic energy and 
chromomagnetic moment. 

The HQET Lagrangian has been studied to order l/niQ (|l8|). For most pur- 
poses, it is sufficient to truncate the theory at £i; higher-order terms are typically 
used to estimate theoretical uncertainties. 

The relation between the meson mass and the 6-quark pole mass, nih, follows 
from Equations |7| and ^, 

run = m, + A-^l±^ + ... (10) 

- Ai — A2 

ruB* = mb + A \- ... . 

Irrih 
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Here Ai and A2 are matrix elements of the quark kinetic energy and chromomag- 
netic moment operators, 

Ai = {B{v)\K{iDfK\B{v))/2mB, 



Ao = (B(v 



B{v))/6mB, (11) 



and A corresponds to the mass of the hght degrees of freedom in the meson (in 
quark model terms, the mass of the light "constituent quark"). The chromo- 
magnetic moment operator is the leading term in the effective Lagrangian that 
breaks the spin symmetry of the heavy quark, and so its matrix element may be 
determined from the B*-B mass splitting, 

A2(mb) ~ 0.12GeV2. (12) 

A and Ai cannot be determined simply by meson mass measurements, although 
on dimensional grounds they are expected to be of order Aqcd and ^qqy), re- 
spectively. The extraction of these two parameters is equivalent to determining 
the pole mass, m^, to relative order 1/m^, and is discussed in detail in Section 
1211. 



2.1.1 OPE and Inclusive Quantities 

The typical momentum transfer in 6-quark decay corresponds to a distance r ~ 
l/m^, whereas the complicated physics of hadronization only arises at the much 
longer distance scale r ~ 1/Aqcd- Because quarks and gluons hadronize with 
unit probability, one would therefore expect the i?-meson lifetime to be the same 
as that of a free b quark in the nib 00 limit. Similarly, any inclusive decay 
distribution (in which all final-state hadrons are summed over) should be given 
by the corresponding free-quark distribution. Similar arguments also apply to 
the cross section for e'^e~ — > hadrons (|l9|) and the r hadronic width (pO|). 

Like the simplifications of HQET, the above argument depends on the sep- 
aration of scales, mf, Aqcdj and one can construct an effective field theory 
description to simplify the analysis. Technically, instead of integrating out heavy 
degrees of freedom as in an EFT, one integrates out the highly virtual degrees of 
freedom via an operator product expansion (OPE) ( ^|2^ , but the anal- 
ysis proceeds in much the same way. Radiative corrections to the heavy-quark 
limit are proportional to powers of as{mb) and are determined by perturbative 
matching conditions onto the OPE, whereas power corrections that scale like 
{Aqcd / mf,)"' are determined by matrix elements of higher-dimension operators. 

For concreteness, we consider the the b ^ u semileptonic width, although 
similar arguments hold for any inclusive decay distribution. From the optical 
theorem, the partial width dT for B decay may be written as the imaginary part 
of the time-ordered product of the weak b ^ u current Ji^ and its adjoint, 

dT ~ Y.^B\Jl^\X){X\J^\B)L^, 

X 

~ l^{B\T{Ji^,j!^)\B)L^,, (13) 
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where = U7^(l — 75)6, L^^, is the lepton tensor, and the summation is over 
all possible final states X. Although the T-product is nonlocal at the scale m^, 
it appears local at the scale Aqcd, and may therefore be written via the OPE as 
a sum of local operators arranged in powers of l/m^: 

TiJii,J^) ~ coif,)KK + "-^KiiDfK + ^^Ka'^'Go^pK + • • • , (14) 

where the coefficients Cj(/i) are calculable in perturbation theory. This is illus- 
trated schematically in Figure |l|. Inclusive quantities may therefore be written 
as a double expansion in powers of AqcD/^i-fe (power corrections, from matrix 
elements of local operators) and as{nrn,) (perturbative matching corrections). 
The OPEs for all inclusive distributions share three important features: 

1. The leading term in the \/rrn, expansion reproduces the free-quark result, 
as argued on physical grounds. In particular, the kinematics are given by 
quark, not hadron, masses. This makes inclusive processes sensitive probes 
of m^. 

2. The 0(l/mfe) corrections to the parton-model result vanish, owing to the 
absence of a dimension-four operator that cannot be removed by the equa- 
tions of motion (^,p2|). 

3. The leading power corrections to the parton result are proportional to the 
matrix elements of {B\b{iD)'^h\B) and {B\bigaf^i,G'^'^b\B) . These are exactly 
the parameters Ai and A2 that entered the relation between nih and in 



Equation IC 



This approach also incorporates an implicit assumption of "quark-hadron du- 
ality," the replacement of the sum over hadron states in Equation ^ with the 
corresponding sum over parton (quark and gluon) states. The validity of this 
assumption clearly depends on how inclusive the observable is. For example, the 
semileptonic 6 — > u decay rate to hadronic states with invariant mass mx < rn.,^ 
obviously vanishes, whereas the corresponding decay rate to free quarks and glu- 
ons is nonzero. However, in regions of such severely restricted phase space, the 
OPE breaks down — the coefficients of both the perturbative and nonperturba- 
tive corrections become large — and there is no question of applying it here. The 
deeper question is whether there are effects due to quark-hadron duality violation 
that do not show up in the OPE, limiting its accuracy even for more inclusive 



distributions (|25, 26, 



Technically, duality violation arises because the OPE is only properly defined in 
the deep Euclidean region, where the particles that are integrated out are highly 
off-shell, whereas for inclusive decays the OPE is performed in the Minkowskian 
region, in which they can go on-shell. Analyticity arguments suggest that extrap- 
olation to the Minkowskian region is valid (^,^^l|), but they do not provide a 
rigorous justification. Reference ( p6D argues that duality violation arises because 
the expansion in powers of AQCD/^ife in the OPE is asymptotic, and duality- 
violating effects do not arise at any finite order in h.QQ,Y:,/mi). Reference ( |25| ) 
argues, based on quark models, that such violations could be large, whereas two- 



dimensional models of QCD (27) suggest that duality violation is a small effect 
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for B decays. Reference (^8|) concludes from a variety of considerations that 
duality violation in semileptonic decays is a negligible effect. 

Because none of these arguments is truly rigorous, quark-hadron duality should 
be tested experimentally. The agreement between the r hadronic width ( pO| ) and 
the OPE suggests that such corrections may not be large, but consistency among 
several different predictions of the OPE would bolster one's confidence that such 
effects may be neglected. 



2.2 NRQCD 

Bound states of a heavy quark and antiquark are nonrelativistic systems in the 
heavy-quark limit. The effective theory to describe such systems is known as 
nonrelativistic QCD (NRQCD). 

Because both NRQCD and HQET correspond to expanding about the niQ 
oo limit, the operators in the two theories are the same. But because the physics 
of a nonrelativistic bound state is very different from that of a single heavy quark 
interacting with light degrees of freedom, the power counting of operators in 
NRQCD differs from that in HQET. 

There are only two important scales in HQET (if we neglect the light-quark 
masses): the heavy-quark mass, mq, and Aqcd- Hence, HQET operators may 
be classified by their order in Aqcd/"iq- In NRQCD, the dynamics depends on 
two additional scales: the heavy-quark momentum, pq = nigv, and its kinetic 
energy, Eq = ^rnQv"^ (where v is the relative three-velocity of the two heavy 
quarks). Because Eq/mQ is the same order in the nonrelativistic expansion as 
Pq/rriq, the relevant expansion parameter in NRQCD is not 1/mq but rather 
the heavy-quark velocity v. The velocity scaling rules of NRQCD operators have 
been worked out (|29|) . 

A formulation of NRQCD was proposed by Bodwin, Braaten, and Lepage 
(BBL) (|30|), and the analogous theory for electromagnetism, NRQED, had been 
developed earlier by Caswell & Lepage (|3l|) . The BBL form of the NRQCD 
Lagrangian is 

^heavy = -^0 + SC, (15) 

where the leading term is 

V'(x) is a two-component Pauli spinor, and there is a corresponding term for the 
antiquark field x{x). The leading corrections are suppressed by O(f^) relative to 
Cq and are given by 

6C = ^i;\■D^fi; + ^^|;\■D.gB-gB.T>)^|; (17) 
smg smq 

+ -^V^(«D X gE-gEx iB) • + -^tp^gB • a)^ + .... 
8m^ 2mq 

The important difference between Equations 15-17 and Equations |^-^ is that, in 
HQET, the heavy-quark kinetic energy is treated as a perturbation, whereas in 
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NRQCD it is leading-order. The correct treatment of the kinetic term is impor- 
tant for describing nonrelativistic Coulomb exchange. For bound states or states 
near threshold, v ~ Os, and Coulomb-exchange graphs proportional to (os/w)™ 
must be summed to all orders. At leading order in v, this is equivalent to solving 
the nonrelativistic Schrodinger equation. Because NRQCD incorporates relativis- 
tic corrections through higher-dimensional operators, it provides an elegant tool 
for calculating the relativistic corrections to nonrelativistic quantum mechanics, 
including a correct treatment of radiative corrections and renormalization (which 
are much more subtle in any other formulation, such as Bethe-Salpeter). 

In several proposed alternative formalisms for NRQCD, the scales m, mv, and 
mv'^ are explicitly separated, which is convenient for many purposes 



32, 



2.3 Lattice QCD 

In lattice field theory, the spacetime continuum is replaced by a discrete lattice 



(see Reference (34) for reviews of lattice QCD). This introduces an ultraviolet 
cutoff, p = n/a (where a is the lattice spacing). The lattice QCD Lagrangian 
contains, for example, discrete differences, which replace the derivatives of con- 
tinuum QCD. Because discretization errors are short-distance effects, we can use 
field-theoretic methods to separate them from the long-distance QCD dynamics. 
As first proposed by Symanzik (|3^), EFT can be used to study the effects of the 
lattice discretization. 

Symanzik showed that the lattice theory can be described by a local effective 
Lagrangian. The leading-order term is simply the continuum QCD Lagrangian. 
However, the effective theory also contains higher-dimension operators (which are 
accompanied by the appropriate power of the lattice spacing) : 

£eff = >Cqcd + aCi + a^C2 + . . . , (18) 

where i2QCD is the (Euclidean) continuum QCD Lagrangian, 

Cqc^ = -q{If) + m)q . (19) 

The Ci {i = 1,2,.. .) contain operators of dimension d = 4 + i, and all discretiza- 
tion effects of the lattice theory can be described by the Ci. The coefficients 
of the operators in the Ci depend on the underlying lattice theory. Hence, the 
discretization effects of the lattice theory are organized as a power expansion in 
the lattice spacing, a.Q The operators typically scale with the momenta of the 
participating particles in the process, which is Aqcd for light quarks (neglect- 
ing quark masses). Hence, we need aAgcD ^ 1 (or more generally, ap <C 1) 
for the expansion of Equation |l^ to be well-behaved. For example, the leading 
discretization effects of the Wilson action (|3^) are described by 



Ci = ciqa^^F^^'q (20) 



^The situation is a bit more complicated because of the implicit a dependence of the coeffi- 



cients through their dependence on the QCD parameters, as and niq. 
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in the effective tlieory. A priori, several operators contribute at 0{a). However, 
for the matching between the lattice and effective theories, we need to consider 
only on-shell quantities (37, 38). We can therefore use the equations of motion 



to reduce the number of operators that contribute at any given order of o. This 
leaves only one dimension-five operator in the Ci term for the Wilson action. In 
summary, the leading lattice-spacing artifacts of the Wilson action are O(aAQCD) 
(neglecting quark-mass effects), and we see that the a ^ limit recovers contin- 
uum QCD from lattice QCD. 

This formalism naturally leads us to consider improved formulations of lattice 
QCD, where operators are added to the lattice QCD Lagrangian so that the 
coefficients of the leading-order corrections in £efr vanish. This procedure is 
known as Symanzik improvement. For example, we can add a discretized version 
of the dimension-five operator in Equation ^ to the Wilson action to obtain a 
lattice action [the Sheikholeslami-Wohlert action (^)] that is correct to O(a^), 
and hence leaves Ci = 0. If n-loop perturbation theory is used to match the 
lattice and effective theories, then the 0{a) improved action will only be correct 
up to terms of order a"+^a. 

The Symanzik formalism implicitly assumes that arug <C 1, which can be seen 
explicitly by considering, as an example, the energy-momentum relation obtained 
from the Wilson or Sheikholeslami-Wohlert actions (Hol): 



where 



E^ip)=ml + ^p^ + 0{p^), (21) 

7712 



^ = l + 0(aV) (22) 
7n2 

and mi and m2 (the rest and kinetic masses, respectively) are functions of the 
lattice bare mass, mo. Although this lattice artifact is suppressed by a^, it is 
large when am ~ 0(1). The lattice artifacts that lead to the breakdown can be 
identified as higher-order terms in the effective Lagrangian of the form (■jqDq)"', 
with n > 2. These terms can be eliminated by the field equations (^,^2[): 

i^oDor ^{Tng+^f-Dr. (23) 

We see that when am^ ~ 0(1); these higher-order terms in the effective La- 
grangian are no longer small, and the expansion of Equation [l^ breaks down. 

With currently practical lattice spacings, am;, > 1. If we don't want to wait 
until we have the computational power to reduce the lattice spacing to am^ <C 1, 
we must modify the above prescription to treat 6-quark mass effects. Because 
TUb ^ AqcD) the 6-quark mass introduces an additional short-distance scale into 
the problem, and we should be able to find an effective-theory framework that 
allows us to lump the short-distance effects from both the lattice spacing and the 
heavy-quark mass into the coefficients of the effective theory. In fact, there are 
three solutions to this problem (for a review of lattice methods for heavy quarks, 
see Reference (^3])). 

The first solution discretizes the continuum effective theory, either HQET 
(Equation |^ for the i?-meson system or NRQCD (Equation 15) for the bb sys- 



tem. The static theory (44) simply discretizes the leading term in the HQET 
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Lagrangian, Equation |8|: 

/:static = V'^^o'V , (24) 

where Djf* is a discretization of the continuum covariant derivative. The leading- 
order lattice NRQCD Lagrangian is (E 



/ a(2) \ 

/:lnrqcd = ^o'* - ^ J , (25) 

where A^^) is a discretization of the Laplacian operator. The NRQCD propa- 
gators are determined by an evolution equation, which is computationally much 
simpler to solve than the matrix inversion required by the Dirac propagators. As 
in the continuum, relativistic corrections can be added to the leading-order term 



using a discretized version of the correction operators in Equation |17|. Lattice- 
spacing errors can be corrected in a similar fashion by adding new operators to 
the lattice Lagrangian. In this procedure, similar to Symanzik improvement, the 
coefficients of the correction operators are obtained by matching to the continuum 
theory. NRQCD is a nonrenormalizable theory, as evidenced by the power-law 
divergences of some of the coefficients of the NRQCD operators (^5|). As a result, 
the continuum limit (a — > 0) cannot be taken explicitly. Instead, lattice-spacing 
errors are controlled by adding more terms to the lattice NRQCD Lagrangian 
until these errors are sufficiently small. 

The second solution starts with the observation that the Wilson action has 
the same heavy-quark symmetries as continuum QCD (|40|). Indeed, in the limit 
aniQ — > oo, the Wilson action reduces to the static limit (^4|), which corresponds 
to the leading term in HQET. Hence, instead of matching our relativistic Wilson 
action to continuum QCD, we can match it to continuum HQET (|^,|l|). The 
difference between the matching for B and bb systems is the power counting of the 
operators. In this prescription, the operators of the continuum effective theory 
are the same as in the usual HQET, as defined, for example, in Equation |^, 
albeit with different coefficients. All discretization effects are again contained in 
the coefficients of the operators of the effective Lagrangian. Hence we have a 
modified HQET of the form 



with 
and 



■^HQET — ^0 + C'l + ... , (26) 
C'q = hy{iD ■ V — mi)hy (27) 

4 = ^hv{iDfK + ^Ka'^^G^pK • (28) 

Note that we have incorporated a rest-mass term into the leading-order term of 
the HQET Lagrangian. This term is usually omitted in the standard HQET pre- 
scriptions (see Equation ^). We may add it to our modified HQET Lagrangian, 
since it has no effect on the dynamics and affects the mass spectrum only ad- 
ditively (||,||,||,|l|). The coefficient of the kinet ic term in Ci matches the 



coefficient of the usual HQET if the lattice quark mass is adjusted so that 

m2 = mQ . (29) 
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The kinetic term in CJ^ can easily be matched nonperturbatively, if Equation 29 is 
imposed on hadron masses. The above prescription demonstrates exphcitly that 
the difference between m\ and 7712 is a lattice artifact that has no effect on the 
dynamics of the system. These arguments were recently confirmed in a numeri- 
cal simulation of heavy-quark systems using a relativistic 0(a)-improved lattice 
action (|47|). In order to adjust the coefficient of the chromomagnetic interaction, 
c^, to its continuum counterpart, we need the 0(a) improved Sheikholeslami- 
Wohlert lattice action. At present, the chromomagnetic operator is matched at 
tree level only. 

In this framework, lattice artifacts arise from the mismatch between the coef- 
ficients of higher-dimensional operators of the two effective theories. Just as in 
the Symanzik formalism, lattice artifacts can either be reduced by brute force 
(taking a — > 0) or by adding higher-dimensional operators to the lattice action. 
However, when considering higher-dimensional operators for building improved 
actions, we must now allow for "nonrelativistic" operators in our lattice action, 
where spacelike and timelike operators have different coefficients (pO|). 

The third solution is based on the observation that it is possible to match 
relativistic lattice actions with am ~ 0(1) to Equation 15 at the cost of adding 
an additional parameter to the lattice Lagrangian ( 40 , 42 ) . This parameter sep- 
arates timelike and spacelike operators starting at dimension four. With this 
additional parameter, one can impose m\ = and recover the relativistic 
energy-momentum relation. The disadvantage of this method is that it adds 
an additional parameter that must be adjusted (either perturbatively or non- 
perturbatively) to recover Equation |l^. This method was recently tested in a 
numerical simulation ( p7D with good results. 

All three solutions discussed above yield lattice QCD formulations that treat 
heavy-quark-mass effects correctly and allow a systematic analysis of both dis- 
cretization and heavy-quark-mass effects. Most numerical simulations of heavy- 
quark systems are based on either the first or the second method discussed above. 



3 QUARK-MASS DEFINITIONS 

Like any parameter in a Lagrangian, the quark mass is a renormalized quantity 
and must be appropriately defined by some renormalization condition. In princi- 
ple, any renormalization condition is allowed; however, some are more convenient 
to use in a particular calculation than others. In this section, we discuss several 
renormalization schemes for the 6-quark mass and highlight their advantages and 
disadvantages for specific problems. 

3.1 The Pole Mass and Renormalons 

The simplest quark-mass definition is the pole mass, defined as the solution to 

^ - m - S(p,m)| 2=„2 =0, (30) 
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where m) is the self energy, in terms of which the fuh quark propagator is 

iS{p,m) = - — -. (31) 

p — m — 1j[p, m) 

For simpHcity, we denote the 6-quark pole mass by mf, in this paper. 

The pole mass is the simplest definition to use in HQET and NRQCD, and 
is related to the meson mass m b via the 1 /mf, expansion in Equation [l^. It is 
gauge-invariant ( p8| ) and infrared- finite (49,5^). However, despite its simplicity, 



the pole mass has the disadvantage that the perturbation series relating it to 
physical quantities (such as the decay width) is typically very poorly behaved. 
For example, the relation between and the semileptonic 6 — > u width is given 
at two loops by (51, 52) 



Tib^XJ,,) = m^(,_2Al^^ (32) 



+ (3.39 - 3.22/?o) ^ ) +••• 

12^5 



as(mf,)^ 



1927r3 ^ ' ' 



where 



/3o = 11 - fnj (33) 

is the coefficient of the QCD beta function,^ Uf = 4 is the number of light 
flavors, the dots denote terms of higher order in or l/nib, and we have taken 
as{mij) = 0.22. The large two- loop term indicates that the perturbation series is 
poorly behaved, even though as(mf,)/7r <^ 1. 

This can be seen in a different way by using the prescription of Brodsky, Lep- 
age & Mackenzie (BLM) (|53|), in which the 0{alPo) term is used to determine 
the appropriate scale for as(//) in the one-loop term. In the BLM prescription, 
the O(a^/3o) piece of the two- loop correction (arising from vacuum polarization 
graphs) is absorbed into the one-loop term by a change of renormalization scale; 
the resulting scale is taken to be the appropriate renormalization scale of the 
process. In the case of Equation |3^, this leads to a scale 

/ 2(3 22) \ 

/XBLM = m^exp f " ~ 300 MeV. (34) 

Since this scale is much lower than the typical momentum transfer in the problem, 
it indicates a deeper problem. 

Bigi et al. and Beneke & Braun have demonstrated (^,55) that this sickness 



persists to all orders in perturbation theory. These authors showed that any per- 
turbation series that relates a physical quantity to the pole mass has an intrinsic 
ambiguity of relative order Aqcb /n^b because perturbation theory is only asymp- 
totically convergent (|56[). Because can be determined only by its relation to 
a physical quantity, this translates into an ambiguity of the same order in the 
^Some authors define (3o with an additional factor of 1/4. 
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definition of tlie pole mass itself. The poor behavior of the series Equation 32 
even at second order appears to be a reflection of this. 

This type of intrinsic ambiguity in perturbation theory is known as an infrared 
renormalon (^7\j. Physically, it arises from the low-momentum region of loop 
integrals where QCD is strongly coupled. Infrared renormalons are ubiquitous in 
QCD perturbation theory and do not signal an inherent limitation of the theory; 
in a consistent OPE, there is always a nonperturbative matrix element that enters 
at the same order in Aqcd/Q as the ambiguity in perturbation theory (where Q 
is the hard momentum scale in the process), and the sum of perturbative and 
nonperturbative effects is well-defined (57). 

In contrast, the leading nonperturbative term in the expression for the semilep- 
tonic 6 — > u width enters at 0(AQ(^£,/m^). Hence, it cannot absorb the ambiguity 
in the series in Equation which means that there is an inherent ambiguity of 
O(Aqcd) in the pole mass. The pole mass is particularly sensitive to infrared 
physics because it is defined as a property of an unphysical on-shell quark. The 
better-defined mass parameters we discuss in the next section are renormalized 
at momentum scales much greater than Aqcd- Hence, they are insensitive to 
long-distance physics and do not have this problem. 

Rather than discuss the formal theory of infrared renormalons (see Refer- 
ence (|57[)), we illustrate their effect with an explicit example. Because it is 
not feasible to calculate to arbitrarily high order, little has been established rig- 
orously about the asymptotic behavior of QCD perturbation theory. However, 
a qualitative understanding may be obtained by considering the class of terms 
proportional to o'^'^^Pq (the higher-loop analogues of the BLM term). This se- 
ries is usually simple to compute, since it corresponds to replacing the gluon 
propagator at one loop with the geometric series shown in Figure |2| and then 
substituting nf ^ — 3/3o/2. Note that this is not a well-defined expansion: There 
is no /3o — *■ oo limit of QCD, and there is no reason that these terms should dom- 
inate perturbation theory. However, the series provides a tool to examine high 
orders of perturbation theory, and barring any miraculous cancellations, we ex- 
pect the conclusions we draw from this subset of graphs to remain valid in QCD. 
The techniques to perform the bubble sum were presented in Reference (^) . The 
series of Equation |3^ continues (retaining only terms of order a"/3Q~^) as 

Gl\Vu,?ml 



1927r3 

3 



1 -2.41^+ 40.5 (35) 



-820 (^i^) + 19192 (^i^) - 533495 (^i^)' 
+1.75 X 10^ (^iM) ' _ 6.71 X 10« (^^) ' 



7 

+ 



= ^' ' ^ (1 - 0.15 - 0.11 - 0.09 - 0.09 - 0.11 - 0.15 - 0.25 - . . . . 

1927r'^ 

After several terms, perturbation theory begins to diverge, as expected. The 
ambiguity in the series is the same size as the smallest term in the series, which 
a more formal analysis shows to be of order AgcD/^T-fe- 

A similar situation arises for the relation between m;, and any other physical 
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quantity. For example, the first moment of the hadronic invariant mass spectrum 
in semileptonic b ^ u decay is related to the pole mass by (p 



= 0.202^^+ 3.151 + 51.91 (^^)' (36) 



m 



+940.52 ( ^f^y + 19347.5 ( + ^ A + . . . 

7 A 

= 0.014 + 0.015 + 0.018 + 0.023 + 0.032 + + . . . 

10 niB 

which shows no signs of converging. The dependence of Equation ^ on is 
implicit in the A term. The two-loop term corresponds to a BLM scale /^blm — 
0.03 ~ 100 MeV, and once again an infrared renormalon ambiguity is present 
at 0{Aqcd /rrib). 

However, eliminating rrih (or equivalently. A) and expressing the two physical 
quantities F and {sh) in terms of one another results in the much better-behaved 
relation 

M = -0.135^^ - 0.601 f + 1.56 f (37) 

mi, TT \ TT J \ TT J 



+148.1 ('^iM V + 4923. ( + + . . . 



TT J \ TT J 10 TUB 

10 TUB 



~ -0.0086 - 0.0024 + 0.0004 + 0.0026 + 0.0051 + — + 



where A^' = m^ — (1927r^F/G|.|14fcp) 5 is a physical quantity of order Aqcd- [Falk 
et al. ( |59| ) used this to define the "decay mass" rri^ = uib — A^' + 0(Aqcd/''7t.^)i 
which does not suffer from a renormalon ambiguity.] The convergence of per- 
turbation theory has improved dramatically. It can easilily be shown that the 
0(Aqcd/"^6) ambiguity has vanished, and the leading renormalon in Equation ^ 
is now at 0(AQ(-,j-,/m^), reflecting the presence of additional unphysical parame- 
ters such as A^ in the OPE for {sh)- The corresponding BLM scale is now 

//BLM = 0.37mf, ~ 1.8 GeV, (38) 

which is significantly greater than before, and is the natural scale for the process 
(recall that the energy must be divided among several final states, so the 
typical momentum transfer in the decay is somewhat less than m^). The poor 
behavior of perturbation theory is therefore an artifact of using the pole mass as 
an unphysical intermediate quantity. 

While this is a problem with the pole mass in principle, in practice there is 
nothing wrong with using it as an intermediate quantity, as long as it is used 
consistently. However, the presence of the renormalon ambiguity in mj, results 
in pole mass determinations that strongly depend on the order of perturbation 
theory used in the calculation. Hence, a practical disadvantage of this approach 
is the difficulty of estimating the theoretical uncertainty in mf,. 
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3.2 The MS Mass 

In general, a mass parameter renormalized at a scale /U is insensitive to physics 
at longer distance scales. A short-distance mass is not plagued by the infrared 
problems that afflict the pole mass. An example of such a mass definition in 
a momentum-subtraction renormalization scheme is the Georgi-Politzer mass, 
defined at the spacelike subtraction point —p^ = m? (^). However, most analytic 
perturbative calculations are performed using dimensional regularization. As a 
result, the most common short-distance mass definition is the MS mass mb(;u), 
which is defined by regulating QCD in dimensional regularization and subtracting 
the divergences in the modified minimal subtraction scheme. The renormalization 
scale /i of a short-distance mass is typically chosen to be of the same order as 
the characteristic energy scale Q of a process, since perturbation theory typically 
contains terms proportional to q;^(/x) log"^{Q/fi), which are otherwise large. The 
renormalization group may be used to relate a mass renormalized at one scale to 
that at another scale, summing all the logs of this form to all orders. 

The relation between the MS mass and the pole mass is known to 0{a^) (|6l| , 
11,11): 



1 + |^ + (1.562/?o- 3.739) f^) (39) 



mh 3 vr V ^ 

+ (l.4686^^ + 0.2548/?i + 0.3905/3o - 29.94) (—^ + 0(af), 



where 



and we use the notation 
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(3i = 102 - —Uf (40) 



mb = mbiiTib), (41) 

where rribijnb) is the MS mass renormalized at the scale ^ = fnb{p,) and as = 
<^s{^b)- All light quarks have been treated as massless in this expression. The 
complete rric dependence of the O(a^) term is known (pl|), and the O(a^mc) 
terms have been calculated (|6^). 

The MS mass rnb{n) arises naturally in high-energy processes, such as Z ^ bb, 
in which the b quarks are produced relativistically. For example, the contribution 
of the axial current to the decay Z — > bbX is ( |65| ) 

^ ^^aUi + — -6^ i+- (i_2in:::f ) +...} (42) 



+ ... 



In the pole mass scheme, the large logarithm of m^/rn^ makes perturbation 
theory artificially badly behaved; choosing the MS mass renormalized at a scale 
^ ~ mz eliminates the large logarithm. 

However, the MS mass is less useful for heavy quarks at nonrelativistic energies. 
Trading the pole mass for nib only slightly improves the apparent convergence of 
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perturbation theory for the semileptonic b width at two loops (|5ll), 

Tib^XJu) = roM^ (l + 4.25^+ 4.58 (^)%o + ...) 

~ ToM^ (1 + 0.30 + 0.19 + ...), (43) 

corresponding to a BLM scale /iBLM — 0.12 m;, ~ 500 MeV. However, the 
asymptotic behavior of the series improves, and indeed it can be shown that 
the O(Aqcd) ambiguity vanishes in this relation. 

Because the typical momentum transfer in semileptonic b decays is somewhat 
less than nih, it might be argued that the appropriate mass to use is mh{^), where 
/i < 772-6 . Indeed, Reference (^) stressed that the appropriate renormalization 
point for b ^ u semileptonic decay is mi^/n, where n = 5 is the power of in the 
total semileptonic width. However, as emphasized in Reference (|66|), the MS mass 
is not a useful quantity when renormalized at a scale ^ < m^. From the effective- 
theory perspective, this is easy to see; the MS mass is defined in full QCD, treating 
the b quark as fully dynamical. This is the appropriate theory to calculate the 
running of the MS mass from some high scale down to mf,. At the scale = m^, 
however, the effective theory changes from QCD to HQET. It therefore makes no 
sense to lower // beyond rrn, in full QCD, and in HQET the pole mass does not 
run. Renormalizing fnh{fj,) below this scale simply introduces spurious logarithms 
that have no physical significance, and therefore do not improve the convergence 
of perturbation theory. Thus, although the MS mass is at least well-defined, it is 
not a particularly useful quantity to relate to low-energy observables. 

3.3 Threshold Masses 

Because of the shortcomings of the pole and MS masses for describing the physics 
of nonrelativistic heavy quarks, several alternative mass definitions have been sug- 
gested, which we group here under the term "threshold masses" (following Ref- 
erence (|67|)). These are mass definitions for which the ©(AqcD/^ife) renormalon 
is absent, but which have better-behaved perturbative relations to properties of 
nonrelativistic heavy quarks than the MS mass. One may define an arbitrary 
number of sensible mass parameters, but several definitions have become popular 
in the literature, and we now consider them in turn. 

3.3.1 The Kinetic Mass 

The kinetic mass w-^'"(^j) introduced by Bigi and collaborators 
defined by introducing an explicit factorization scale ^f, and subtracting the 
physics at scales below fij from the quark-mass definition. 

More explicitly, the kinetic mass is defined by considering various sum rules 
for semileptonic b ^ c decay in the small velocity (SV) limit (^), mb,mc ^ 
mb — rric 3> Aqcd- In this limit, the charmed meson is produced with vanishingly 
small recoil {rric ^ — m-c), but there is still enough energy transfer to produce 
a large number of excited charmed states {rrib — rric ^ Aqcd)) so that the decay 



68, 6e, 
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may be treated inclusively. In this limit, one can derive sum rules that relate A 
and Ai (and higher-order terms if required) to weighted integrals of the spectral 
function. By putting an explicit cutoff /i/ on the integrals, Bigi et al. define 
a cutoff A(/xj) and Xi{pf) that determine the kinetic mass via (compare with 
Equation pX| ) 



ruB = mf'^ifif) + A{^if) 



Ai(^/) 



+ 



(44) 



where the A2 term and terms of higher order in l/m^ have been neglected. In 
the limit fj,f ^ 0, the pole mass is regained, whereas for Aqcd <^ fJ-f ^ mb, 
this definition removes the dangerous low-momentum region from the definition 
of m^'°(///) and therefore eliminates the infrared renormalon ambiguity while 
leaving the heavy-quark expansion valid. 

The relation between the kinetic mass and the MS mass is known to two loops 



mf, I 9 ^ 9 




where 



K 



2 \ 6 48 / 144 18 V 2 



X, = log?^-|, X, 
rrih 6 



2/x^ 13 
log ^ - — 
nib D 



(45) 

(46) 
(47) 



and terms of order fi^j/ml and higher have been neglected. 



3.3.2 The Potential-Subtracted Mass 

The potential-subtracted (PS) mass proposed in Reference (|7l| ) has similar prop- 
erties to the kinetic mass, but arises from consideration of the properties of 
nonrelativistic quark-antiquark systems. 

The dynamics of heavy quarkonium are determined by the Schrodinger equa- 
tion, 

(-— + V{r)-E]G{f,0,E)=S^^^r), (48) 
\ rub J 

where E = ^/s — 2mb is the binding energy and V{r) is the static QCD potential. 
This expression includes the total static energy of two heavy quarks at a distance 

r (^,EID> 

^stat(r) = 2m5 + F(r). (49) 

Because this is a physical quantity,^ it is well-defined and should not suffer from 
a renormalon ambiguity. Indeed, the high-order behavior of V{r) has been shown 

^There are also power corrections to V{q), but these are of order Aqpn/g^, so they cannot 
absorb the renormalon ambiguity ([Fsl). 
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to precisely cancel that of the pole mass, so that the combination of Equation ^ 
is well-defined. The infrared sensitivity of the long-distance quark-antiquark po- 
tential exactly cancels that of the pole mass. 

This cancellation is made explicit by eliminating the pole mass in terms of 
the so-called potential-subtracted (PS) mass. The coordinate-space potential is 
defined as the Fourier transform of the momentum-space potential, 



V{r) 



(2^ 



(50) 



As noted in Reference (|7lD, the coordinate-space potential is more sensitive to 
infrared physics than the momentum-space potential because of the contribution 
to the Fourier integral from the region of small \q\, and this region was identified as 
the source of the leading renormalon in V{r). In the PS scheme, this contribution 
is subtracted from the potential and instead included in the mass through the 
definitions 



m 



b^iHf) =mb- 5m{fif), V{r, fif) = V{r) + 2Sm{fif), 



where 



Sm 



1 



V{q). 



(51) 



(52) 



2 y|,-i<^, (2^)3 

Using the PS mass and subtracted potential in the Schrodinger equation thus 
results in a better-behaved perturbation series for the quark mass. 
The relation between m^^ and m;, is known to three loops (|7l|): 



PS 



as(/x)Cir j Q^(^) 



ai - /3o An ^ - 2 



+ 



as(/^) 
47r 



a2 - (2ai/3o + /3i 




+ Pi nn2^-41n^ + 



/4 

^2 



(53) 



where 



ai 
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/4343 ^ 2 

h47r2 

y 162 

'55 

T 



vr 
1 



22 
Y 



- + =fC3)Cl 



1798 



l(iC3]CFTRnf + 
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+ yCs^ CATRUf 



653.71 - 66.354n/ + 1.2346n^ 



(54) 



and Ca = 3, Cp = 4/3, Tr = 1/2. The constant 02 depends on the three-loop 
static potential and was calculated in Reference ([7^). Combining Equations ^ 
and ^ gives the three-loop relation between m^^ and fni,. 

Note that for the appropriate choice l-i-f nih, both the kinetic mass and 
the PS mass may be made to differ from the pole mass by 0{mf,v'^). This is 
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important for power counting in Coulomb systems, since this is the same size as 
the Coulomb binding energy. For both nonrelativistic bb problems and B decays, 
the scale fif ^ 1 GeV gives well-behaved series. 



3.3.3 The 15 Mass 



Both the kinetic and PS masses are defined by introducing an explicit factoriza- 
tion scale /i/ to remove the troublesome infrared physics of the pole mass. In 
constrast, the 15 mass introduced in References (75) and ([7^), which we denote 



here by , achieves a similar goal without introducing a factorization scale. 
However, the renormalon cancellation is subtle in this case, and the 15 mass is 
a well-behaved parameter only if the orders of terms in perturbation theory are 
reinterpreted dTSj). 

The 15 mass is simply defined as one half the energy of the 15 bb state, 
calculated in perturbation theory. To three loops ([77[), 
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+ 



vr 



+ ... 



where 



In 



(56) 



and the other parameters are defined in Equation Note that is renormalization- 
group- invariant . 



In the large /3o limit, the series of Equation 55 contains terms of order 

al,alPo,aif3o,... , 
whereas perturbation theory typically contains terms of order 

This apparent mismatch makes it unclear how the use of the 15 mass will improve 
the convergence of perturbation theory. However, as shown in Reference (^5|), at 
high orders in perturbation theory the coefficient of a"~''^/3o in the continuation 
of Equation |5^ contains terms of the form + . . .-|-1), which exponentiate 

at large n to exp(£) = ^/ {mbasCp)- This factor corrects the mismatch between 
powers of a<j and [3q (at least at large orders in perturbation theory), allowing the 
renormalon ambiguities to cancel between Equation ^ and other perturbation 
series. This observation led Hoang et al. (^) to propose the so-called Upsilon 
expansion. In this approach, terms in the perturbative expansion of the 15 mass 
of order a" are formally taken to be of the same order as those of order a"~^ 
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in other series. To make this manifest, a power-counting parameter e = 1 is 
introduced, and terms of 0(a") in the IS mass are multiphed by e"~^, while 
those in other series are multiphed by e". When combining series, terms of the 
same order in e are combined. 

Using this approach, the IS" mass has been shown to have remarkably well- 
behaved perturbative relations to other physical quantities. For example, the 
6 — > u semileptonic width is given by 

r(6 ^ XJue) = ^^'^g^y""^' (l - - 0.0356|lm + •••)' (^7) 

where we have taken as{mi,) = 0.22, and egLjy; denotes only the terms of order 
enhanced by Pq. 



3.4 Lattice Quark Masses 

Almost all determinations of the 6-quark mass from lattice QCD use lattice per- 
turbation theory to calculate the b quark's self energy. A few general comments 



on lattice perturbation theory are therefore given in Section 3.4.1, which also 
defines the conventions used in the following discussion. 

Before we discuss the specific strategies for 6-quark-mass determinations, we 
briefly review how (light) quark masses are determined from lattice QCD in gen- 
eral. The quark masses are adjustable parameters in the lattice QCD Lagrangian. 
One then calculates a suitable hadron mass on the lattice and adjusts the lattice 
quark mass until the lattice result agrees with the experimentally measured value. 
A suitable hadron is one that is easily and reliably calculable from lattice QCD. 
Because they have fewer valence quarks, mesons are simpler to simulate than 
baryons. Mesons that are stable under the strong interactions are less affected 
by sea-quark effects, the incomplete inclusion of which often causes the dominant 
error. Hence, in light-quark systems, the pion and kaon masses are generally used 
for quark-mass determinations. For the b quark, the most suitable hadrons are 
the (spin-averaged) B and Bg mesons and the T system. 

The procedure described above yields a nonperturbative determination of the 
lattice quark mass as it appears in the lattice QCD Lagrangian, which is used 
in the numerical simulation. Lattice quark masses, though well-defined and free 
of renormalon ambiguities, are not very useful parameters in continuum calcula- 
tions. For light quarks, the relation between the lattice and MS mass is known 
at one loop (for all light-quark actions used in numerical simulations): 



fn{fx) = Z^^{fi)mo (58) 

with 

Z^(/.) = 1 + ^(c(°) +/ln(/xa)2 + ...). (59) 

The value of c^^^ depends on the lattice action, and 7*^ is the anomalous dimen- 
sion at leading order. Several groups have introduced procedures for nonper- 
turbatively determining renormalized quark masses. These procedures avoid the 
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perturbative uncertainty of Equation 58 but must still deal with lattice-spacing 



and other systematic uncertainties. One example is the renormalization-group- 
invariant (RGI) mass as defined in Reference (|78|). One needs perturbation theory 
to obtain m from the RGI mass. Because the relation is known at four-loop order 
and the matching can be done at a very high-energy scale, the conversion of the 
RGI mass to m introduces only a small additional uncertainty. 

There are several different strategies for lattice determinations of the 6-quark 
mass. The first method is similar in spirit to light-quark mass determinations. 
One calculates the kinetic mass^Q the coefficient of the kinetic term in the energy- 
momentum relation. Nonrelativistically, 

E{p)=E{Q) + -^p' + .... (60) 



The dispersion relation of Equation 60 is applied to the hadron system, and the 
lattice quark mass, mo, is tuned until the lattice calculation of Mkin agrees with 
its corresponding experimental result. 

As discussed above, this procedure yields a nonperturbative determination of 
the lattice quark mass, mo, which can be related to the MS mass in perturbation 
theory. This relation is generally known at one-loop order: 

m,{^i) = mo (l + c(i)(^)^^ + ...), (61) 

where c*^^'*(/x) depends, as usual, on the lattice action, and the renormalization 
scheme and scale of the coupling. Numerically, S-^^ is 0(1). 

In the second method, the 6-quark pole mass is determined from lattice QCD 
calculations of the binding energy. For the i?-meson system, we have 

nib = M^""^ - £ , (62) 

where M^^^ is the spin average of the experimentally measured B and B* masses, 
and the binding energy £ is obtained from 

£ = Slat - Eo . (63) 

-Eiat is the binding energy in the i?-meson system calculated from a numerical 
lattice NRQCD or HQET simulation. Eq is the b quark's nonrelativistic self en- 
ergy,0 which depends on the underlying lattice NRQCD (or HQET) action. It 
is a short-distance quantity and hence calculable in perturbation theory. In cal- 
culations with relativistic heavy-quark lattice actions (which contain an explicit 
rest mass term) , Equation |6^ is modified to 

£ = Mi,iat - mi , (64) 

where Mi i^t is the spin-averaged rest mass of the i?-meson as calculated from 



lattice QCD and nii is the lattice heavy quark's rest mass (defined in Section 2.3 ) 
which is again calculable in perturbation theory. 

*Note that this mass differs from the kinetic mass defined in Section 3.3.1. 
In HQET language, this term is also known as the "residual mass, 
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Both £^0 and mi have been calculated to one-loop order in perturbation theory 
(8^,81), for example, 

Eq = ^ap + 0{al>) . (65) 

Calculations of the two-loop corrections for both Eq and rrii are currently in 
progress (H. Trottier, private communication), and we expect these results to 
become available soon. Once the two-loop results are available, it should be 
possible to estimate the three-loop correction, using a numerical technique that 
has been succcessfully used to determine the three-loop correction for the static 
self energy ( p2[ ) (see Equation ^ below). The coefficient eo in Equation ^ 
depends on the lattice NRQCD action used in the numerical calculation, and on 
the renormalization scale (and scheme) of the coupling. It also depends mildly 
on the lattice quark mass. The power-law divergences present in both £^0 and the 
numerically calculated Eis,t cancel, leaving the binding energy £ divergence free. 

In the static limit, the b quark's self energy was until recently known only to 
two-loop order (|83|, [79|). The three-loop coefficient has now been calculated by 
two groups ( |8^ , |8^ ) using different numerical techniques: 



171 00 

-^0 



1.0701 ap{OM/a) + 0.117a|, + (3.56 ± 0.50)a^ + 0{aj,) . 



(66) 



The Hf dependence of the two- loop coefficient is known; reported here are the 
values at = (in the quenched approximation). The coefficient of the Op 
term is only known at nj = 0. 



Comparing Equation 62 to Equation 10, we see that a lattice calculation of the 
binding energy, £, can be used to determine the HQET parameters. A, Ai 
A2) (HD: 

•^1 , / 2\ 



(and 



£ = A 



2mki, 



+ 0{l/m' 



(67) 



Here l/2mkin is a short-distance coefficient defined from the heavy quark's energy 



momentum relation as in (for example) Equation 21 or Equation 6C. The term 
A2 drops out of this equation, since £ in Equation ^ is calculated from the spin 
average of the B and B* masses. It can be determined separately by considering 
the B*-B mass difference. 

For the bb system. Equations 62-^ must be modified to account for the two 
heavy quarks in the bound state: 



rrib 



i(M^"P 



and 



£ = Slat -2E0OT £ = M 



l.lat 



2mi . 



(68) 



(69) 



In either case, the pole mass determined from Equation 62 (or Equation |6q) can 
be converted to any other continuum mass defined in the previous subsections, 
using the relations given there. The renormalon ambiguity present in the pole 
mass manifests itself in the perturbative expansion of Eq (or mi). If one uses. 



for example. Equation 39 to relate the pole mass to the MS mass, then the 



renormalon ambiguities in Equation 39 and 62 will cancel, leaving a well-defined 



perturbative expansion. In order to make this cancellation explicit, one should 
use the same coupling in both equations. 



26 



3.4.1 Lattice Perturbation Theory 

The previous section shows that lattice perturbation theory is often a necessary 
ingredient in quark-mass determinations. However, perturbative expansions ex- 
pressed in terms of the bare lattice coupling are not well-behaved. For example, 
the static self energy of Equation ^ expressed in terms of the bare lattice cou- 
pling, aiat, takes the form 

= 2.1173 aiat + 11.152 afat + (86.2 ± 0.5) al, + 0{ai,) . (70) 



The perturbative series in Equation |70| looks very divergent, with increasing co- 
efficients. However, as explained in Reference (|8^), this is an artifact of using 
a poor expansion parameter. The reliability of lattice perturbation theory is 
easily tested by comparing short-distance quantities calculated in Monte Carlo 
simulations to their perturbative expressions. As discussed in Reference (^6|), 
even two-loop predictions fail miserably at reproducing Monte Carlo results for 
short-distance quantities if the perturbative results are expressed in terms of the 
bare lattice coupling. If instead the perturbative results are expressed in terms 
of a renormalized coupling (such as ay or aj;jg), then perturbative and Monte 
Carlo results are in good agreement. The accuracy of perturbative predictions is 
further improved if the scale at which the coupling is evaluated in the perturba- 
tive expansion corresponds to the typical momentum scale of the gluons in the 
physical quantity. 

We follow the procedure suggested in Reference (p6|), which is particularly 
convenient for lattice perturbation theory and has been shown to produce reliable 
perturbative estimates. We define the coupling ap from the expectation value of 
the plaquette (the smallest Wilson loop on the lattice): 

-ln(Trf/p) = — ap(3.40/a)(l - 1.1909ap). (71) 

The coupling ap is defined so that the perturbative expansion of the plaquette in 
terms of ap contains no higher-order terms. The definition of ap is designed to 
coincide with the coupling defined from the heavy-quark potential in momentum 
space, 

through one-loop order, 

ap = av + 0{al.) . (73) 
The scale q* at which the coupling is evaluated is defined as 

where f{q) is the integrand of the quantity that is evaluated at one loop. This 
definition for setting the scale is very similar in spirit to the BLM procedure for 
continuum perturbation theory. 

This procedure turns the perturbative expansion of Equation ^ into the much 
better-behaved Equation 
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4 DETERMINATIONS OF THE 6-QUARK MASS 



4.1 The bb System 

The T system has historically been an important source of information on mi,. As 
discussed in Section |l|, potential models provide a good fit to the observed spec- 
trum of bb resonances and give model-dependent determinations of m,;,. More re- 
cently, there have been model-independent determinations of m;, from the masses 
of the low-lying resonances as well as from the near-threshold behavior of the 
e'^e~ — > bbX cross section. All of these determinations use the fact that bb states 
are nonrelativistic, and therefore can, to first order, be described by the non- 
relativistic Schrodinger equation. Effective field theory is then typically used to 
calculate both relativistic and nonperturbative corrections to this limit. 

In this section, we discuss determinations of from the T(1<S') mass using 
both calculations based on perturbation theory and those from lattice QCD. We 
also discuss determinations of nih from T sum rules. 



4.1.1 The T(IS') Mass 

In the heavy-quark limit, the bb pair form a nonrelativistic Coulomb bound state 
with dynamics determined by the Schrodinger equation. Thus, in the heavy-quark 
limit it is straightforward to determine nih from the spectrum of T mesons. For 
sufficiently heavy quarkonium, 

mr(ns) = ^r(ls) + ^^rinS) ' C^^) 

where ^71^(^5) calculable in perturbation theory and '5Afx(^5) denotes the non- 
perturbative contribution to the T mass. This immediately gives the result 

ml^ = 4.73 GeV - I^M^P^^ (76) 

or 

frib = 4.21 GeV it nonperturbative corrections. (77) 

The trick is to determine the size of the nonperturbative corrections. 

Potential-model studies indicate that the T system is far from a Coulomb 
bound state; the radius of even the lowest lying states sits squarely in the confining 
potential, which suggests that nonperturbative effects are important for these 
states. [More recently, Brambilla et al. ( ^7|) showed that this may merely reflect 
the use of the poorly defined pole mass in the potential model. These authors 
find that perturbation theory expressed in terms of a short-distance mass gives a 
good fit to the T(IS'), ^(l-P), and T(2S') levels of bottomium.] Voloshin ( ^8|) and 
Leutwyler (|89|) showed years ago that the leading nonperturbative corrections to 
the heavy-quark limit cannot be absorbed in a nonperturbative potential. Rather, 
in the heavy-quark limit a bb bound state is a compact object of size <^ I/Aqcd, 
and nonperturbative effects correspond to the interaction of the QCD vacuum 
with the multipole moments of the bound state. Furthermore, the correlation 
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time for the background gluon field ~ 1/Aqcd must also be much greater than 
the dynamical timescale for the quarkonium state ~ l/rrihv'^, where v ~ as{mbv) 
is the typical velocity of the b quark in the meson. 

In this limit, the nonperturbative correction SM^^^g--^ may be written as a series 
in terms of local condensates (88, 89|,pC|): 



oo 9/a 9\2/. 9 \ 2r+2 

sM?u.^a^[^) . PS) 

The obvious difficulty with treating the bb system in this approach is that the 
inequality 

mbalinibv) > Aqcd (79) 

does not hold; both sides are of order a few hundred MeV. Nevertheless, because 
the leading operator only arises at 0(Aq^£,), taking the OPE at face value gives 
a rather small correction to the T(15') mass from the first condensate, 

^Mrns) = ^^rrib ^^^c^^^^A MeV, (80) 

using {<^\asG''^'''Gl^^ = 0.05 GeV^ (0) and evaluating at the Bohr radius of 
the 15 state. This is a reasonably small contribution, even though the inequality 



of Equation 79 is not well satisfied. The more realistic scaling 

nibV > rubv^ ~ Aqcd (81) 

is much more difficult to describe theoretically. The corresponding corrections to 
the heavy-quark limit are given by condensates that are nonlocal in time, about 
which very little is known (p^). 



Because Equation pQ is not expected to dominate the series in Equation 78 for 



physical values of the 6-quark mass, it is not clear how sensible this estimate of 



the nonperturbative corrections is. References (75,93) and (p4| ) take Equation 8C 
as indicative of the size of nonperturbative corrections and treat it as an estimate 
of the theoretical error rather than a correction. Beneke & Signer (p^) then 
determine 

mP(2 GeV) = 4.58 ± 0.08 GeV, (82) 

which corresponds to 

Mfo = 4.24 ±0.09 GeV. (83) 

Pineda ( |9^ ) estimates contributions of 0(Aq(jj-,) and 0(Aqqj-,) condensates and 
finds that their effects largely cancel (underscoring the poor behavior of the ex- 
pansion), giving 

m6 = 4.21 ±0.09 GeV. (84) 
Brambilla et al. (^7|), including the effects of the charm-quark mass but ignoring 



nonperturbative corrections entirely, find 



mft = 4.19 ±0.03 GeV. 



(85) 
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4.1.2 Lattice Calculation of the bb spectrum 

At present, determinations of the 6-quark mass from the bb spectrum using lattice 
QCD are limited by our knowledge of the heavy-quark self energy at finite quark 
mass (away from the static limit). The 6-quark mass is known only at one-loop 
order in perturbation theory, which leaves an uncertainty of ~ 100 MeV in m;, 



The NRQCD group (^, ^7j) has calculated the 6-quark mass by using a 
lattice NRQCD action for the b quarks that is correct through 0{v^) and O(a^). 
They use gauge configurations in the quenched approximation at several lattice 



spacings (98), as well as gauge configurations with Uf = 2 light sea quarks using 
two different actions for the light quarks: staggered (|9^ ) and 0{a) improved 
Sheikholeslami ( |97|) fermions. They find 



mb{mb,nf = 0) = (4.28 ± 0.03 ± 0.03 ± 0.10) GeV (86) 
mb(mb,nf = 2) = (4.26 ±0.04 ±0.03 ±0.10) GeV , (87) 

where the first error is statistical, the second error includes an estimate of higher- 
order relativistic and discretization effects, and the third error estimates the 
perturbative uncertainty. From their study of the lattice-spacing dependence, the 
authors find that nib changes within 30 MeV — in agreement with their systematic 
error estimate. The difference between the quenched and Uf = 2 result is smaller 
than their statistical errors. They have also studied the dependence of on 
the sea-quark mass (|97|), with somewhat heavy sea quarks. The change of mf, is 
negligible compared with the statistical errors when the sea-quark mass is varied 
between 2m ^ and m,. 



4.1.3 Spectral Moments 



Moments of e~^e — > bbX distribution were proposed to determine heavy-quark 
masses from quarkonium more than 20 years ago ( 99| , |100| , 101), but there has been 
a resurgence of interest recently ( p^[T03t[l0^ , |l05| , p^ , [T07| ,|9^,[l0^,|6l 
particularly with the advent of NRQCD technology, which simplifies the calcula- 
tion of subleading corrections for large moments. The hope is that because one 
integrates over several resonances, the sum rules are less sensitive to nonpertur- 
bative corrections than are the properties of the individual bound states. 
Consider the correlator of two electromagnetic currents of bottom quarks. 



where 

j^(x) ^ Kxh^bix). (89) 

n(g^) has poles in the complex plane at the location of bb bound states, and a 
cut on the positive real axis corresponding to the continuum. By analyticity, the 
n'th derivative of n(g^) is therefore related to an integral along the cut 
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while the optical theorem relates the imaginary part of the vacuum polarization 
n(g^) to the total cross section for e+e^ ^ 7* ^ 66 + X, 



Rb{s) = UnQllmUiq'^ = s + ie) (91) 



where 



CTpt "^35 3 

The resulting sum rule relates derivatives of the vacuum polarization Il{q^) to 
the experimentally measurable moments of the total cross section for e~^e~ — > 66 
pairs: 



n! d{q 
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On dimensional grounds, the left-hand side is proportional to so a precise 

value of m,}, may be determined for large values of n. 

In practice, n must be neither too large nor too small. Because the experimental 
measurement of R{s) is very poor in the continuum, n must be large enough that 
the moment is dominated by the first few T resonances, whose properties are 
known quite accurately [with current experimental data, a ±50 MeV experimental 
error on mf, requires n > 6 (|93[)]. On the other hand, as n increases, the sum rule 
is dominated by low-momentum states near threshold. Hence, nonperturbative 
effects become increasingly important in the calculation of the left-hand side of 
Equation |9^ as n increases. These may be determined by expanding the product 
of currents in Equation ^ in an OPE. This gives a series of the same form as 
Equation in which the m^w and mbV^ are replaced by the typical momentum 
and energy of the 66 states that dominate the moment. As before, the leading 
nonperturbative contribution comes from the gluon condensate {agG""^^ G^j^) . 
This was estimated to be a < 1% effect in m;, for n < 20 (100, 102). Because 



the relative size of this term grows like n^, nonperturbative corrections to the 
moments are negligible for smaller values of n. 

This approach may underestimate the size of nonperturbative corrections. The 
characteristic energy and momenta for the nth moment are (|105|, P^) 



E ^ rrih/n, p~2mf,/-v/re. (94) 

The characteristic energy becomes O(Aqcd) for n ~ 10, which suggests that 
nonperturbative effects should be larger than this simple estimate. Hoang ( |105[ ) 
therefore argues that a reliable extraction of rrib cannot be obtained from moments 
n > 10. 

With this caveat in mind, for appropriate values of n the left-hand side of 



Equation 93 is perturbatively calculable. However, for moments large enough 
for the continuum to be neglected, the intermediate 66 states are produced near 
threshold, and perturbation theory breaks down just as it does for the calcula- 
tion of bound-state properties. The Coulomb-enhanced terms, which are pro- 
portional to as/v for bound states, are now proportional to as\/n ~ 1 and 
must be summed to all orders. This is done in the same way as for 66 bound 
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states. Solving the nonrelativistic Schrodinger equation sums all terms of the 
form {as^/n)^, while perturbative corrections (proportional to powers of as) and 
relativistic corrections (proportional to powers of 1/y/n) are parametrically the 
same size and are calculated using NRQCD technology. For details, see Refer- 
ences ( |10^ , [10^ , |108| , |6^, |107|, p^) . 

The current state of the art in these calculations is NNLO (next-to-next-to- 
leading order, or relative order 1/n). Several groups have performed the analysis 
(see Table These results are all consistent with one another, with errors esti- 
mated at the ±100 MeV to ±30 MeV level on the threshold masses, and somewhat 
larger errors on nib. 

The different groups use different approaches to determine their central values 
and error estimates (Reference ( p3|) compares the methods), but in no case does 
the nonrelativistic expansion seem particularly well-behaved. References ( |105| ) 
and ( |111| ) fit for the pole mass mi,, so the result has large variation at different 
orders because of the infrared renormalon discussed in Section 3^ . Perturbation 
theory is greatly improved when written in terms of the short-distance kinetic, 
PS, or IS" masses ([T06|,P,|107|,[9^), but the NNLO corrections are still disturbingly 
large, as illustrated in Figure ^ (from Reference (93)). As shown in the figure, 
the NNLO corrections shift m^^ by about 250 MeV compared with the NLO 
result for a renormalization scale ~ 2 GeV. This is about the same size as 
the shift in nib from LO to NLO, so perturbation theory does not appear to be 
converging well. Beneke & Signer (93) trace this behavior to the corresponding 
poor convergence of perturbation theory for the leptonic width of the T(15'), 
which dominates the sum rules for large n. 

Given this poor behavior, it is difficult to estimate the theoretical error in nib 
reliably. A conservative approach would be to vary the renormalization scale 
between mb/2^/n ~ 1.5 GeV and 2mb/\/n ~ 6 GeV and determine the error from 
the difference between the NLO and NNLO calculations. As Figure ^ shows, 
this corresponds to a variation in m^^ from ~ 4.38 to ~ 4.84 GeV, or about a 
±230 MeV uncertainty. Beneke &: Signer (|9^ ) argue that, because the calculation 
becomes unstable below /i ~ 2 GeV, it is more appropriate to determine the 
uncertainty from the variation in the NNLO result alone when fi is varied over 
the range 2 GeV < fi < 6 GeV, which yields their quoted scale uncertainty of 
±100 MeV. 

Melnikov &i Yelkhovsky ( |107 ) and Hoang (105,106,34) find similar behavior for 
perturbation theory for individual moments but have different ways of reducing 
the theoretical error. Melnikov & Yelkhovsky (107) interpret the first two terms in 
the nonrelativistic expansion as the first two terms of an alternating series. This 
can be converted via an Euler transformation to a more convergent series, from 
which they extract a smaller uncertainty. [However, these authors use larger 
moments, n = 14-18, whose reliability has been questioned ( [105 , 93).] It is 
difficult to prove that this rearrangement of the perturbative series is valid beyond 
the order at which it is currently calculated. In the absence of a more rigorous 
theoretical argument, it is not clear that the corresponding theoretical uncertainty 
is reduced by this procedure. In another approach, Hoang ( |105| , [106| , |6^ notes 
that the behavior of the nonrelativistic expansion is much improved when one 
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performs simultaneous fits to multiple moments from n = 4 to n = 10. This 
improves the stability of the fit greatly; both the shift in m{, from NLO to NNLO 
and the renormalization-scale dependence of the result are much smaller than for 
single moments, which accounts for the smaller theoretical uncertainty quoted in 
these results. However, the reason for this improved behavior is not immediately 
clear, and the resulting theoretical error estimate depends on the correlations 
between the different experimental uncertainties in the masses and widths of the 
low-lying resonances. Once again, it is difficult to prove that this approach really 
gives a better-behaved perturbation series. 

A poorly converging perturbative expansion frequently indicates that a para- 
metrically large class of terms has not been resummed. The situation is similar 
to the calculation of the production cross section a{e~^e~ ttX) near thresh- 
old, where the analogous calculation has been performed to NNLO (for a recent 



review, see Reference (112)). As with the bb sum rules, the NNLO correction to 
the cross section is roughly the same size as the NLO correction, and there is a 
large renormalization scale uncertainty in the NNLO result. It was shown for ti 
production that perturbation theory is much improved by the use of the renor- 



malization group in NRQCD to sum logarithms of the form log™ v (|33| , 113). 
For bb sum rules, the analogous renormalization-group-equation (RGE) calcula- 
tion would sum terms of order a™ log™ n; it would be interesting to see if this 
gives a similar improvement. 

In the absence of such a calculation, the poor behavior of perturbation theory 
suggests that the theoretical situation is not as stable as the values in Table 
H suggest. Given these issues, it is probably prudent to assign a conservative 
theoretical error to these determinations, at least until the poor convergence of 
perturbation theory is better understood. 

Finally, Kiihn & Steinhauser (109) avoid the issue of resumming Coulomb cor- 
rections by considering low moments (n < 4) for which fixed order perturbation 
theory is appropriate. As noted earlier, for such low moments the extracted 
value of rrib is sensitive to R{s) in the continuum regime where it is poorly mea- 
sured. The authors replace the experimental measurement of R{s) above the 
resonance region {^/s > 11.2 GeV) with its QCD prediction. This greatly reduces 
the uncertainty of nib since the perturbatively estimated uncertainty on the QCD 
prediction of R{s) close to threshold is significantly smaller than the experimental 
uncertainty. The authors find 

mfe(mb) = 4.21 ± 0.05 GeV (95) 

which is consistent with the determinations from higher moments. However, 
unlike the determinations from higher moments, the uncertainty in this result 
depends strongly on the assumption that the QCD prediction for R{s) (with 
perturbatively estimated errors) is valid very close to threshold, where it has not 
been well measured. 
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4- 2 The B System 

Unlike the bh system, the B system does not become perturbative in the heavy- 
quark hmit, since the size of the hadron is still determined by nonperturbative 
physics. Hence, lattice QCD is best suited for determinations of from the 

5- meson spectrum. On the other hand, as discussed in Section |2|, inclusive quan- 
tities such as the S-meson semileptonic width are dominated by short-distance 
(r ~ l/'Tift) physics, and are therefore sensitive to and perturbatively cal- 
culable in the heavy quark limit. Power corrections to the heavy-quark limit 
that scale like {^qicd / mh)^ may be parameterized in the framework of HQET, 
allowing a precision determination of rrn,. 

Below we discuss determinations of rrib from the i?-meson spectrum and from 
inclusive B decays. Typically, results in the B system are quoted for A and Ai 
rather than a better-defined threshold mass. Because A is only defined order by 
order in perturbation theory, it is important to work consistently in Us- In this 
section, we therefore distinguish A extracted at one and two loops. 



4.2.1 Lattice Calculations of the i?-meson Spectrum 



The best determination of the 6-quark mass from lattice QCD calculations of the 
S-meson spectrum comes from calculations of the binding energy in the static 
limit, because the static self energy is known at three-loop order. 

At present, there are two independent determinations of the 6-quark mass from 
the static binding energy. The first one uses the results from numerical calcula- 
tions of the S-meson system using static h quarks ( 114 , 115| ). Allton et al. (114) 
use an 0{a) (tree-level) improved action for the light valence quarks. They per- 
form their calculations at several lattice spacings and use the quenched approx- 
imation {uf = 0). Gimenez et al. ( |115| ) obtain their results from numerical 
simulations with Uf = 2 light Wilson quarks. Based on the numerical results of 
References (114) and (|115|), the authors of References (|79|,|116|), and (115) obtain: 



mb{mb,nf = 0) 
rnb{mb,nf = 2) 



(4.30 ± 0.05 ± 0.05) GeV 
(4.26 ± 0.06 ± 0.07) GeV . 



(96) 
(97) 



The first error in Equation ^ is dominated by the uncertainty in the lattice 
spacing but also contains the statistical error of the numerical simulation and 
the effect of the uncertainty in a^. The second error in Equation 96 is the 
perturbative error in both Equation and Equation 66. Similarly, the first 



error in Equation 97 is due to statistical and other systematic lattice errors, 
while the second error is the perturbative error. Because the nj dependence of 
the three- loop term in Equation |6^ is unknown, at present, this error is larger 
than in the quenched case. 

The second determination (117) uses the numerical results of References (|118|) 



and (119), which present numerical calculations of the i?-meson spectrum using 
lattice NRQCD for the b quark and an 0{a) improved action for the valence light 
quarks. The results of Reference (118) are obtained in the quenched approxima- 
tion, whereas the results of Reference (|119|) come from numerical simulations with 
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Uf = 2 staggered sea quarks. After extrapolating the results of Reference ( |118| 
to the static limit, Collins (|117|) obtains a 6-quark mass of 



rnb{mb,nf = 0) = (4.34 ± 0.04 ± 0.05) GeV . (98) 

The first error combines the statistical uncertainty with the uncertainty in the 
lattice spacing, and an estimate of residual discretization effects of O^o^Aqqj^). 
The second error is the perturbative uncertainty. Collins ( |117| ) estimates that 
the inclusion of sea quarks lowers rrib by about 70MeV, having used two-loop 
perturbation theory to compare rnb{rnb,nf = 2) with rnb{rnb,nf = 0). 

Both determinations use preliminary results for the coefficient of Op in Equa- 



tion 66, which have a larger uncertainty than the final result shown in Equa- 



tion I63. Hence, the perturbative uncertainties given in Equation |96| and Equa- 
tion ^ are slightly overestimated. In both determinations, a mild residual lattice- 
spacing dependence was observed, which is roughly equal to the statistical uncer- 
tainty, 10-30 MeV. Hence, in order to resolve the residual lattice-spacing depen- 
dence, the statistical accuracy of the numerical simulations must improve. This 
should be feasible with currently available computational resources. The error 
due to using the static limit is 0(AQCD/"^^,), and affects rub at the ~ 1% level. 
This error will be removed when two- and three-loop results for the heavy-quark 
self energy at finite quark mass become available. The nj dependence is similar 
in size to the perturbative uncertainty. Numerical simulations with Uf = 3 sea 
quarks are needed to bring this error completely under control. 

Heitger & Sommer suggest a new method for determining the 6-quark mass 
based on the nonpertubative methods developed by the ALPHA collaboration (|120|). 



The authors obtain a preliminary result for mf,(4GeV) in the quenched approxi- 



mation, which corresponds to mb{mb) = 4.48 ± 0.13 GeV (121). Error analysis is 
in progress. 



As discussed in Section 3.4, the determination of the binding energy £ away 



from the static limit can be used to extract both A and Ai by fitting the bind- 



ing energy to Equation |67|. Kronfeld & Simone (85) used results of numerical 
simulations of the i3-meson system (as obtained in References ( |122 , 12^ , 124 )) to 
determine these parameters. They find 

A = 0.68l°:?2 GeV (99) 

and 

Ai = -(0.45 ±0.12) GeV^ (100) 

where the error includes statistical and lattice-spacing errors. These results come 
from numerical simulations that were performed in the quenched approximation, 
for which no solid error analysis exists at present. The results are based on one- 
loop perturbation theory, in which the coupling is evaluated at the scale q* ~ 1/a, 
determined from Equation ^ (which is similar to the BLM scale). The errors 



in Equations 99 -IOC do not include an estimate of the perturbative uncertainty. 
Earlier determinations of A and Ai were based on numerical simulations of lattice 
HQET ( |T25D . 
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4.2.2 QCD Sum Rules 

The first QCD-based determinations of mf, from tlie B system relied on the QCD 
sum rules of Shifman, Vainshtein & Zakharov (SVZ) (|9ll). In this approach, 
correlators of hadronic currents are studied at intermediate distances (of order 
~ 1 GeV~^). At this energy scale, one can either expand the product of currents 
in an OPE as a sum of local operators, or evaluate the currents by inserting 
complete sets of states between the operators. For sufficiently small separation, 
the OPE is dominated by the lowest-dimension operators; for sufficiently large 
separation, the sum over intermediate states is dominated by the lowest-lying 
hadronic states. QCD sum rules rely on a "window" of distance scales for which 
both approaches hold (in practice, the convergence is improved by the use of Borel 
transforms, and the window is in the Borel parameter), and use local duality to 
model the contributions of all but the lowest-lying states with the perturbative 
QCD prediction. In this way, local duality is used to relate hadronic matrix 
elements (from the sum over states) to a small number of condensates (from the 
OPE). 

QCD sum rules have been remarkably successful at describing many low-energy 
features of QCD and have been used extensively to study properties of heavy 
quarks, including rrif, (126, 127, 125, 128] ) (for reviews of recent results, see Ref- 



erence (|130| )). Because of their well-defined heavy-quark limit, QCD sum rules 
may be used to determine the parameters of HQET. However, unlike the other 
approaches discussed in this article, they are not strictly model-independent; it 
is difficult to quantify the systematic errors inherent in the assumption of lo- 
cal duality in the intermediate window. As a result, the sum rules cannot be 
consistently improved by calculating higher-order corrections. 

With this caveat in mind, the HQET sum rule determination of A gives (125, 

400 - 600 MeV Reference (128) 
A(lloop) = <( (101) 
570 ± 70 MeV Reference (OT 



The sum rule for Ai is more problematic. Two calculations report confiicting 
results: Ai = -0.5±0.2GeV2 ( |l3l| ) and Ai = -0.1 ± 0.05 GeV^ (|l32|). The 



discrepancy has been traced to the the contribution of an off-diagonal matrix 
element involving an excited pseudoscalar B meson (|69D, but it is currently un- 
resolved. 

Several extractions do not explicitly use the heavy-quark limit. From the sum 



rules for the pseudoscalar two-point function, Narison (127) obtains the two- loop 
result 

rub = 4.59 ± 0.06 GeV ^mb = 4.05 ± 0.06 GeV, (102) 
where the errors do not address the assumptions inherent in the sum rules. 



4.2.3 Inclusive Moments 

As discussed in Section ^, inclusive differential decay rates of heavy hadrons may 
be used to determine nib without the additional assumptions required by QCD 
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sum rules. The first attempts to extract HQET parameters from inclusive decays 
combined the rates for B — > Xce(/i)Pe(^) and D — > Xu^s^^'e. decay (133, 134) as well 



as i? ^ XcTVr ( |135| ), but these extractions rely on the validity of the heavy-quark 
expansion for inclusive charm decays. As a result, they have large errors. 

With the advent of precision measurements in recent years, direct determina- 
tions of A and Ai may be made purely in the B system. Of the infinite number of 
inclusive quantities, we focus on three that have been widely used to determine 
mf,: the shape of the charged-lepton energy spectrum in i? — > Xcii^i, moments of 
the hadronic invariant mass spectrum in — > Xcii'i, and moments of the photon 
spectrum in B ^ Xgj. Power corrections to these quantities are all known to 
0{l/m^), and the one-loop perturbative corrections are also known. In addition, 
the 0(q^/3o) (BLM) piece of the two-loop term, which is expected to dominate 
the two-loop result, is known for most of these observables. Since there are more 
free parameters at 0(l/m^) than current data can constrain, estimates of the 
sizes of these terms are typically used to estimate the theoretical uncertainty in 
the extraction of A and Ai from a given observable. 

For consistency in comparing different approaches, the estimating technique 
proposed in Reference (|136| ) is often used. In this approach, the parameters at 
0(l/m^) are independently varied over a range of ±(500 MeV)^ and the region 
containing 68% of the points in the A, Ai plane is taken as an indication of 
the theoretical error. This approach allows the relative theoretical errors from 
different observables to be compared in a consistent way, but the arbitrariness of 
this procedure should be kept in mind when interpreting the experimental results. 
Changing the range of variation of the parameters to, for example, ±600 MeV 
would increase the theoretical error by ~ 70%. Figure ^ compares the relative 
theoretical uncertainties in the (A,Ai) plane obtained by this technique from 
hypothetical measurements of moments of the charged-lepton energy spectrum, 
hadronic invariant mass, and photon energy spectrum in i? ^ ^sl, discussed in 
the following sections. 

In addition to the experimental determinations of A and Ai, it is possible to 
derive purely theoretical constraints. By using sum rules that relate exclusive 
decay form factors to Ai and A2 and demanding that the exclusive rate be less 
than the inclusive rate, Bigi et al. (^) derived the constraint 

- Ai > 3A2 ~ 0.36 GeV^ ± . . . , (103) 



where the ellipses denote radiative corrections. The authors (38) use a defi' 



nition of Ai(^) in which the matrix element in Equation ^ has been defined 



with a cutoff as discussed in Section 3.3.1. This differs from the definition 



of Al in dimensional regularization, and the radiative corrections to the bound 



of Equation 103 differ in the two schemes. Equation 103 in the MS scheme is 



considerably weakened by radiative corrections (139); the best bound quoted is 
— Al > 0.01 GeV^, and the two-loop corrections were large enough to challenge 
the validity of perturbation theory. The radiative corrections to the bound on 
Ai(^) for ^ ~ 1 GeV are expected to be better behaved (|69D. 

The difference between the two schemes is formally higher-order than the ex- 
pressions we consider in this section. Within experimental errors, the bound in 
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Equation |103| appears to be satisfied. 

The Charged-Lepton Energy Spectrum: Moments of the charged-lepton 
energy spectrum of the form 

( \ f^f dV 



Mr = I ^^dEe"^^' 

were suggested ( |140| ) as a sensitive probe of the mass difference mf, — rric. How- 
ever, because of the large background from secondary leptons (for example, from 
charm decays), it is difficult to measure the charged-lepton spectrum over the 
full kinematic range. [CLEO (|141|) imposed a cut Eg > 1.4 GeV, and modeled 
the spectrum at lower energies.] To avoid the model dependence of the low-E'^ 



region, Gremm et al. ( 142 ) proposed the following observables: 



^ dr ,^ f^r-' dr ,^ 

Ei—dEi / TErdEe 

n — 1.5 GeV dEj „ _ Jl.TGcV dEj nnc;\ 

'^^-dE, ' ^dE, 



1.5 GeV dEi 7l.5GeV dE£ 

Gremm and colleagues calculated expressions for i?i and R2 to 0(l/m^) and 
0{as) ( |142| ), and calculated the 0(l/m^) corrections ( |136| ) and the 0{alPo) 
terms ( |143| ). A correction due to the boost of the leptons must also be taken into 
account (|142| ). Ri and R2 constrain similar linear combinations of A and Ai, 
so the corresponding solid ellipse in Figure ^ constrains one linear combination 
much better than A or Ai individually. 

The Hadronic Invariant Mass Spectrum: Moments of the form 

1 /""^B dr 

{{SH - mD^r ) = - {SH- mD^r-—dsH, (106) 
r Jo dsH 

where sh is the invariant mass of the final-state hadrons and rfiD = (m/3+3m^)/4 
is the spin-averaged meson mass, have been proposed (|59|). At the parton level, 
the invariant mass sh of the final hadronic state in semileptonic b ^ c decay is 
fixed at , so positive moments of sh — rnj^ only get contributions in the OPE 
at 0{as) and 0{l/mh), making them a good probe of the power corrections in the 
OPE. Falk et al. calculated expressions for the first two moments, {{sH—rnj-,)) and 
{{sh — ^i))^) (59), to 0(l/m^) and a^Po, while Gremm & Kapustin calculated 



the 0(l/mg) terms (|T36D 



Falk et al. ( [59| ) used the experimentally measured branching fraction to the ex- 
cited states Di and D2 to put lower bounds on the first two moments, excluding 
a region of the A — Ai plane. The early DELPHI result of 34 it 7% semileptonic 
branching fraction to these states gives a strong bound of A > 410 MeV, but 
later measurements that put this branching fraction closer to 11% weaken this 
bound considerably, to A > 100 MeV. More useful is the direct measurement of 
these moments by the CLEO collaboration, made possible by reconstructing the 
missing neutrino. This introduces a subtlety because the technique of neutrino 
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reconstruction requires placing a lower cut of 1.5 GeV on the charged-lepton en- 
ergy; Reference ( |137] ) presents the complete theoretical prediction, including the 
effects of the lepton cut. 

The first moment of the invariant mass spectrum constrains roughly the same 
linear combination of A and Ai as Ri and R2, discussed above. The second 
moment constrains a roughly orthogonal linear combination, but the theoretical 
prediction is very uncertain because of the effects of 0(l/m^) terms in the OPE, 
so the dashed ellipse in Figure ^ usefully constrains only the same rough linear 
combination of parameters as Ri and i?2- 



The B Xsj Photon Spectrum: Comparison of the measured weak ra- 
diative B Xg^ decay rate with theory is an important test of the standard 
model. In contrast to the decay rate itself, the shape of the photon spectrum is 
not expected to be sensitive to new physics (|144| , |145| ) but instead serves as an 
additional inclusive quantity that may be used to determine A and Ai. 

The photon spectrum in the large-energy region {E^^ > 1.5 GeV) is dominated 
by the operator 



O7 



167r2 



(107) 



with small contributions from four-fermion operators. For lower photon energies 
there is a large background from nonleptonic decays, so it is necessary to restrict 
measurements to large photon energies. Kapustin & Ligeti ( |144| ) studied moments 
of the form 



M(")(i?o) 



Eo 



^7 j!^^^7 



£^0 



dEj ^ 



(108) 



At tree level, measurements of the first moment and variance of the photon energy 
directly determine A and Ai, respectively: 



m!^'\eo) 

Mf (^o)-(M«(i?o))' 



niB — A 



'12 V 2 / 



1 -|- (5mi 
2 



5m 



/2Eo 
/2Eo 



+ 



(109) 



V nib 



25mi 



2Eo 
nib 



+ 



where 6mi and 5m2 are perturbative corrections, and the ellipses denote terms 
of order a^, a^A^coMb' and A^co/^^fe- 

Corrections to these results of order a^/?o ( |146|) and (AqcD/w-fe)^ (138) have 
been calculated, and the estimated theoretical error is shown as the dot-dashed 
line in Figure ^ An important additional uncertainty in these results arises from 
the effects of the photon cut Eq. Only moments of the photon spectrum are 
calculable in the OPE; the precise shape is determined by the nonperturbative 



parton distribution function of the b quark (147) and so the effects of the cut 
must be modeled. Simple models for the distribution function suggest that a 
photon cut oi Eq ^ 2 GeV is sufficient for the dominant uncertainty in A to 
be set by the terms, whereas the corresponding uncertainty in Ai is very 



model-dependent ( 138| , 145 ). Hence, only A is likely to be determined accurately 
by the photon spectrum. 
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Experimental Results: From CLEO measurements (tabulated in Reference (|148| 



Gremm et al. (142) determined the experimental values of Ri and R2, 

P = 1.7831 GeV, R^""^ = 0.6159, (110) 

which lead to the values 

A(l loop) = 0.39 ± 0.11 GeV, Ai = -0.19 ± 0.10 GeV^, (111) 

where the uncertainty corresponds to the experimental la statistical error. The 
estimated uncertainty due to 0(l/m^) terms is larger than this ( [1361 ). The 
O(a^/3o) terms calculated in Reference ( |143D shift the central values slightly, 

A(2 loop) = 0.33 GeV, Ai = -0.17 GeV^. (112) 

More recently, Ri and R2 were calculated in the IS scheme by Bauer & Trott 
(|l4g| ). Using the values in Equation 109, they find 

Ai5 = 0.47 ±0.10e ±0.08t GeV, Ai = -0.16 ± O.ll^; ± 0.08t GeV^, (113) 

where A15 = niB — ml'^ , corresponding to ml^ = 4.84 it O.IO^; it 0.08^. This 
gives the MS mass 

mb = 4.31 ±0.13 GeV, (114) 

adding the experimental and theoretical errors in quadrature. The CLEO collab- 
oration (150) has recently measured these moments, and a preliminary anal- 



ysis gives, for decays to electrons, Ri = 1.7797 ± 0.0007stat ± 0.0007sys and 
R2 = 0.6173 lb 0.0016stat i 0.0014sys (with compatible results for muons). This 
in turn gives m = 4.81 ± 0.09e (electrons) and ml^ = 4.87 ± O.IIe muons , 
where only the experimental error is quoted. 

In 1998, the CLEO collaboration published preliminary measurements of the 
first and second hadronic invariant mass moments, as well as the first two lepton- 
energy moments M^^^ and M^^^ (|151| ). The values of A and Ai obtained from 
the hadronic invariant mass moments differed from the values obtained from the 
lepton-energy moments, although the discrepancy was only at the ~ 2a level. 
Two possible issues are the use of the second hadronic invariant mass moment, 
which is not theoretically well-controlled ( |137] ), and the model-dependent extrap- 
olation of the data into the region Eg < 0.6 GeV ( |152| ). 

More recently, two new published CLEO measurements of the hadronic invari- 
ant mass and photon energy spectra (153, 154) give 



{sH-rn%) = 0.251 ± 0.023 ± 0.062 GeV^ 
{(sH-mlf) = 0.639 ± 0.056 ± 0.178 GeV^ (115) 

and 

(E^) = 2.346 ± 0.032 ± 0.011 GeV 
{{E^ - {E.,)f) = 0.0226 ± 0.0066 ± 0.0020 GeV^. (116) 

Note that the photon spectrum was measured only down to E^/ = 2.0 GeV, and 
the moments fitted to theoretical predictions that included the photon-energy 



40 



cut ( |146| ) . Because of the theoretical uncertainties discussed above for the second 
moments, robust constraints are obtained only from the two lowest moments 
(Figure |5|) . This leads to the values 

A (two loop) = 0.35 ± 0.07e ± O.IOt GeV 

Ai = -0.236 ±0.071^ ±0.078r GeV^ (117) 

(where the theoretical error includes an estimate of the operators as well 



as the scale dependence of a^), in good agreement with Equation 112. 

As we have discussed, it is difficult to assign a perturbative error to the extrac- 
tion of A because of its inherent ambiguity; it is better to work with a well-defined 



threshold mass.^ For the 15" mass, this gives (see also Reference ( 155D ) 

ml^ = 4.75 ± 0.07e GeV (118) 

corresponding to 

mft = 4.22 ± 0.07s ± 0.05i/„3 GeV (119) 

where the first error is the experimental uncertainty and the second error is an 
estimate of the uncertainty due to 0(l/m-^) corrections obtained from the width 
of the ellipse in Figure 4. This result is in good agreement with the results from 
the T system. The perturbative uncertainty is not included in the error, since it 
would require directly fitting the moments in the 15 scheme to the data, which 
has not yet been done. Nevertheless, the theoretical error appears comparable 
to the theoretical error from the T system, though the sources of theoretical 
uncertainty are completely different. 



4-3 High-Energy Determinations 



The DELPHI and ALEPH collaborations at LEP ( 1156 , 157 ) have analyzed e~^e~ 
annihilation into heavy-quark jets, which depends on rrih. Using the NLO pre- 
dictions ( [L58D , they find 



mh{mz) = 2.67 ± 0.25,tat ± 0.34had ± 0.27thy GeV (DELPHI) (120) 
mbimz) = 3.27 ±0.22,tat±0.22exp±0.38had±0.16thy GeV (ALEPH) 

(where "had" denotes hadronization uncertainties) . Although the precision is not 
comparable to that of low-energy determinations, this first measurement of 
away from threshold provides a statistically significant check of the running of 
mb(//). Both results are in good agreement with the predicted QCD running of 
the low-energy masses determined from the T system. For example, running 
the DELPHI result down to nif, gives 

mbifnb) = 3.91 ± 0.67 GeV, (121) 

in agreement with other determinations (but with significantly larger error). 



^This is not the case for the corresponding error in \Vcb\ quoted in Reference (153), since the 



renormalon ambiguity cancels between physical quantities. 
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5 CONCLUSIONS 

All the precision determinations of m;, discussed in this review are currently 
dominated by theoretical errors. Table |2| summarizes the model- independent 
extractions of rrib from the different approaches discussed in this paper. We 
assign conservative theoretical errors to the results, which encompass the ranges 
given in the literature. 

The determination of mf, from the perturbative calculation of the T(IS') mass 
is limited by our knowledge of the nonperturbative corrections. The associated 
error is difficult to estimate and to improve, since the corresponding OPE does 
not appear to converge well. 

Sum rules for spectral moments of bb production are less sensitive to nonper- 
turbative effects and have the potential to give very accurate determinations of 
nib. The main issue in present determinations is the poorly behaved perturba- 
tive series, which might be improved via a renormalization-group resummation. 
In the absence of such a calculation, the perturbative uncertainty is difficult to 
estimate reliably. 

The TTife determination from inclusive S-meson decays is limited by nonpertur- 
bative effects at ©(Aq^j-j/m^), which are estimated to be at the ±50-70 MeV 
level for the moments we have discussed [although these may be smaller for other 
choices of moments (149)]. It may be possible to better constrain these terms 



from the fits to the different moments, or from lattice calculations. Additional 
uncertainty from violations of quark-hadron duality is also possible. This uncer- 
tainty is difficult to quantify but may be experimentally tested by demanding 
consistency between different moments. Converting the determination of A from 
the photon spectrum in B ^ Xg'j in Equation 116 to a threshold mass will al- 
low the perturbative uncertainty to be accurately estimated, and perturbative 
corrections are expected to be well-behaved. 

With the exception of the static result, determinations oirrib from lattice QCD 
are presently limited by the low order at which the perturbative corrections are 
known. However, perturbative calculations of the two-loop corrections are cur- 
rently in progress, and it may be possible to extract the three-loop corrections 
by using numerical techniques. Hence, we can expect a significant reduction of 
the perturbative uncertainty. In principle, this uncertainty can be reduced even 
further through the use of nonperturbative renormalization prescriptions. At 
present, all lattice calculations suffer from the incomplete inclusion of sea-quark 
effects, which is the most troublesome source of systematic error. However, this 
error appears to have only a small effect on the 6-quark mass determinations. The 
results listed in Table |2| were obtained from numerical simulations with Uf = 2 
sea quarks. Given sufficient computational resources, all systematic errors that 
arise in lattice QCD calculations are in principle controllable and systematically 
reducible. Current progress in numerical simulations of sea-quark effects using 
improved actions ( |159D gives us reason to believe that lattice results with good 
control over all systematic errors may be achieved in the next few years. 

To put these uncertainties into context, consider the precision in rrib presently 
required for CKM physics. Probably the quantity most sensitive to rub is \Vub\, 
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as determined from the inclusive i?-meson decay, b u(.v. This rate is pro- 
portional to m^. Hence, an uncertainty in mf, of ±100 MeV corresponds to a 
±12% uncertainty in the total semileptonic 6 — > u decay width, or a ±6% uncer- 
tainty in \Vuh\- In practice, the sensitivity of \Vub\ to nih may be even stronger 
because of the experimental cuts necessary to reduce the 6 — > c background. For 
example, in a recent proposal for determining \Vub\ through a set of optimized 
kinematic cuts ( 160D , the uncertainty in rrib is the dominant source of theoretical 



error: a ±80 MeV (±2%) uncertainty in results in a ±15% uncertainty in 
\Vub\- Thus, for precision CKM physics, the error on nib should be reduced to 
less than 100 MeV. 

As Table |2| indicates, current determinations are approaching this level, and 
the different results agree very well with one another. Indeed, the variation of 
the central values given in Table 2 is smaller than the uncertainties associated 
with each determination. However, given the discussion above, we feel that it is 
inappropriate to reduce the error in the average. Hence, our best estimate of the 
6-quark mass is 

mbirrib) = 4.24 ± 0.11 GeV, (122) 

where the central value and uncertainty have been chosen to accommodate most 
of the spread in the theoretical determinations. If we run this result up to the 
Z-boson mass, we obtain 

rfib{mz) = 2.94 ± 0.09 GeV . (123) 

in good agreement with the high-energy determinations shown in Equation 117. 
Currently, determinations of from lattice QCD are comparable in accuracy to 
determinations from perturbative QCD, and the prospects for improving on the 
current accuracy of nib determinations from several of these methods over the 
next few years are excellent. 
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Figure 2: The bubble sum, which may be calculated to arbitrarily high order in 
perturbation theory. 



Table 1: Determinations of from high moments of R{s) compared to NRQCD 
at NNLO^ 



n 


Quoted mass (GeV) 


mb{rnb) (GeV) 


Reference 


8... 12^ 


ruh = 4.80 ±0.06 


4.21 ±0.11 


(111) 


14. . . 18^ 


m^'"(lGeV) =4.56 ±0.06 


4.20 ±0.10 


<m 


10 


mP(2 GeV) = 4.60 ±0.11 


4.26 ±0.10 


(HD 


4 . . . 10^= 


ml^ = 4.71 ±0.03 
ml^ = 4.69 ±0.03 


4.20 ±0.06 
4.17 ±0.05 


dH) 

(34) 



^For 



brevity, we have added errors quadratically; see cited references for details on the individual 
sources of uncertainty. In each case, the renormalization scale dependence is the dominant 
source of theoretical uncertainty, and the experimental error is negligible in comfjarison. The 
charm-quark mass nic is taken to be zero in internal loops in all but Reference (p 



References 
'^References (P 



) and ([1 07| ) fit single moments, 
and (p4|) simultaneously fit multiple moments. 
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Figure 3: (2 GeV), extracted from the n = 10 moment in T sum rules at NLO 
(relative order 1/^/n) and NNLO (relative order 1/n) as a function of renormal- 
ization scale fi. (From Reference (p^).) The dark region corresponds to the 
experimental uncertainty on the moment. 



-0.1 
-0.15 

;-o.25 

I 

r -0.3 

-0.35 
-0.4 
-0.45 




0.6 



Figure 4: Relative theoretical uncertainties in A and Ai from (a) the lepton energy 
spectrum in semileptonic B decay (solid line), (6) the hadronic invariant mass 
spectrum in semileptonic B decay (dashed line), and (c) the B Xsj photon 



spectrum (dashed-dotted line), from References (|137] ) and (|138| ). 
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Figure 5: Determination of A and Ai from the hadronic energy spectrum and 
photon spectrum in semileptonic and radiative B decays, from Reference (|153| ). 
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Table 2: Summary of model-independent determinations of m}){mbY 



System 


Method (Section) 


Caveat 


frlbimb) (GeV) 


T 


IS" mass 


nonperturbative terms 


4.23 ±0.11 




sum rules 


poor convergence 


4.20 ±0.10 




lattice QCD 


sea-quark effects 


4.26 ±0.11 


B 


Inclusive moments 

(i) lepton spectrum 

(ii) photon energy/ 
hadron invt. mass 


duality 


4.31 ± 0.13 

4.22 ± 0.09 ± 0{al) 




lattice QCD (static limit) 


sea-quark effects 


4.26 ± 0.09 




OUR DETERMINATION 




4.24 ±0.11 



''The error bars represent our estimate of the uncertainties, based on the discussion in Section 4. 
Theoretical and experimental uncertainties have been added in quadrature. The "Caveat" 
column identifies the most problematic aspect of the theory used in each determination, as 
discussed in the text. The last row gives our best determination ofrnb{rnb), as discussed in the 
text. 



